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EDITORIAL 


La revista Madera y Bosques cumple venturosamente 15 años en este 2009, lo cual es 
para los que contribuimos a su edición, motivo de alegría y acicate para continuar produ- 
ciendo una revista de calidad a la altura de las mejores de su género en el mundo de 
habla hispana. Sin embargo, el entorno nacional y mundial en que este aniversario ocurre 
es bastante complejo y problemático. Los temas de la inseguridad y la violencia azotan a 
importantes estratos de la sociedad mexicana y amenazan con crear un clima de intran- 
quilidad generalizado, a pesar de los mejores esfuerzos de las autoridades. En el ámbito 
económico, el mundo se enfrenta a una crisis de enormes proporciones, como pocas 
veces se han visto en la época moderna, quizá comparable en algunos sentidos con el 
gran ‘crack’ financiero de 1929, y por supuesto, México no es la excepción. En el contexto 
forestal nacional, el sector oficial enfrenta una ‘mini-crisis’ propia, derivada de la incapa- 
cidad confesa de la presente administración de cumplir con las metas de reforestación 
que ella misma se habla marcado como compromiso. Todo esto describe un panorama 
que a simple vista podría parecer catastrófico, pero afortunadamente, aún hay razones 
para el optimismo. Es cierto que la inseguridad derivada del tráfico de drogas asuela a 
importantes sectores sociales, pero afortunadamente, como se reconoce tanto en el 
sector gubernamental como en la propia sociedad, los puntos críticos no abarcan toda la 
geografía del país y en general, la situación no ha salido de control y se espera que el 
gobierno recupere paulatinamente los espacios que ahora ocupa la delincuencia organi- 
zada. En el ámbito económico, los gobiernos de los principales países han adoptado 
medidas que se espera contribuyan a resolver en plazos no muy largos los aspectos más 
acuciantes de la crisis. Tampoco se pueden ignorar los progresos sociales alcanzados en 
muchos países del mundo, incluido México y otros de Latinoamérica, que si bien no son 
del todo suficientes, sí permiten vislumbrar una mejoría con el concurso de los varios 
grupos y sectores sociales y oficiales. En el tema forestal, existen también datos que 
indican que la situación no es insalvable. A pesar de que no se alcanzó la meta de refo- 
restación en términos de hectáreas, sí se logró un avance sustantivo en términos del 
número de árboles plantados, lo que necesariamente debe ayudar al medio ambiente y a 
la producción forestal en el mediano plazo. Este panorama debe ser visto más que como 
un ambiente de derrota, como un reto importante a vencer y debe motivar a los diferentes 
actores sociales y económicos a sumar fuerzas para salir adelante. En la revista, hemos 
adoptado esa visión y esperamos que nuestra contribución en la difusión de la ciencia 
forestal contribuya con su granito de arena a esa tan ansiada recuperación. 

Los temas abordados en este primer número del año abarcan conceptos por demás inte- 
resantes. El primer artículo aborda el tema de la distribución vertical de aves en un 
bosque templado con diferentes niveles de perturbación. En este trabajo se destaca la 
presencia de diferentes especies de aves en dos zonas localizadas en el estado de 
México, una con bosque de pino poco alterado y otra con bosque mixto perturbado. Se 
caracterizó y comparó la distribución de aves en estratos verticales de la vegetación y se 
desarrollaron análisis específicos para describir la estructura vegetal de ambas zonas. La 
metodología descrita aquí permitiría detectar y evaluar zonas perturbadas en bosques 
templados de otras regiones, con la posibilidad de aportar conocimiento que contribuya a 
su mejor manejo y conservación. 

El siguiente artículo se refiere a la aplicación de técnicas multi-criterio en el manejo inte- 
gral forestal. En el estudio se evalúa la capacidad del bosque para producir bienes y serví- 
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cios tomando en cuenta diversas funciones del bosque, así como las interacciones del 
ecosistema. La metodología propuesta permite analizar y resolver un problema forestal 
que involucra el manejo de varios bienes y servicios del bosque, tomando en cuenta la 
opinión de productores y técnicos forestales. Las técnicas multi-criterio de programación 
por compromisos permiten determinar la mejor combinación de las variables que 
respondan al manejo simultáneo de los objetivos buscados y son herramientas que 
pueden ayudar a entender la dinámica de los bosques y el manejo simultáneo de varios 
usos y servicios del bosque en el largo plazo. 

Otro artículo trata de la variación en la cobertura de suelo en un ensayo de procedencias 
de Pinus greggii en otro estado del norte de México. Para conocer la protección al suelo 
que brindan las copas de los árboles en un ensayo de procedencias, se midió supeivi- 
vencia, altura y características de la copa y se calcularon área de proyección de copa, 
área de intercepción lumínica y porcentaje de copa. El estudio muestra evidencias para 
seleccionar las mejores procedencias para proteger el suelo del sitio de plantación. El 
tema de la modelación del volumen fustal de Pinus durangensis en un estado del norte 
de México es tratado en otro manuscrito. Esta especie arbórea es la de mayor uso y distri- 
bución en la región de estudio, sin embargo, el conocimiento para estimar con veracidad 
su volumen es limitado, por lo cual aquí se planteó modelar el ahusamiento de dicha 
especie que posibilite generar una ecuación de volumen comercial que garantice su apli- 
cación operativa. El modelo ajustado arrojó buenos resultados en su capacidad predictiva 
y mediante su integración matemática se obtuvo una ecuación de ahusamiento compa- 
tible con los volúmenes comercial y total. Las ecuaciones generadas son una herramienta 
útil para mejorar la precisión en los inventarios forestales. 

Otro trabajo plantea la modelación matemática del secado convencional de madera. En 
este trabajo se explican los mecanismos de transferencia de masa y calor presentes en 
el secado convencional, se escriben las ecuaciones de transporte de masa y calor, y por 
último se muestran los perfiles de humedad y temperatura para encino europeo Quercus 
pedonculae, sometido a secado convencional. El número se completa con un trabajo 
sobre otro recurso forestal de la mayor importancia, el bambú. Se presenta un método 
alternativo para calcular las propiedades geométricas de las secciones de los culmos de 
bambú que se utilizan en ensayos de flexión para obtener las propiedades mecánicas. Se 
propone un procedimiento utilizando imágenes digitalizadas de las secciones de los 
culmos que se procesan con un programa de CAD efectuando una integración numérica; 
con este método el error se pueden reducir significativamente. 

Esperamos que esta información cumpla con su cometido y que permita a los usuarios 
de estos materiales obtener el mayor provecho de los recursos forestales al aportar ideas 
y metodologías útiles y de aplicación clara. Mejor contribución a la solución de los 
problemas nacionales no se puede hacer en una revista de esta naturaleza. 

Raymundo Dávalos Sotelo 

Editor 
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ARTICULO DE INVESTIGACION 


Distribución vertical de aves 
en un bosque templado 
con diferentes niveles de perturbación 

Vertical distribution of birds in a températe forest 
with different perturbation levels 

Saúl Ugalde-Lezama 1 , Juan Ignacio Valdez-Hernández 2 , 
Gustavo Ramírez- Valverde 3 , José Luis Alcántara-Carbajal 1 

y Juan Velázquez-Mendoza 2 

RESUMEN 


De febrero a junio del 2004, la presencia de diferentes especies de aves fue estudiada en dos 
zonas localizadas en la Estación Forestal Experimental Zoquiapan (EFEZ), Estado de México: una 
con bosque de pino poco alterado (ZOO 1) y otra con bosque mixto perturbado (ZOQ 2). Se empleó 
el método de recuento en puntos con radio fijo para el monitoreo de aves; así como los métodos de 
cuadros empotrados y cuadrantes con punto central para el muestreo de la vegetación. Se caracte- 
rizó y comparó la distribución de aves en estratos verticales de la vegetación, mediante análisis de 
regresión Poisson y por frecuencia de observación; además, se desarrollaron análisis específicos para 
describir la estructura vegetal de ambas zonas. La distribución vertical de las aves siguió cuatro 
patrones: bajo (B), alto (A), alto-bajo (AB) e indeterminado (I). La abundancia de individuos por estrato 
tuvo un valor diferente entre estratos y zonas; en ambas, el estrato supenor tuvo una mayor abun- 
dancia que el inferior y medio. Sin embargo, el número de individuos por especie, registrado en cada 
estrato, fue mayor en el bosque perturbado. El Valor de Importancia Relativa (VIR) de la vegetación 
confirmó que ZOQ 1 es un rodal dominado por Pinus hartwegii y algunas arbustivas de Penstemon 
gentianoides y Senecio asteracae. ZOQ 2 es un bosque mixto conformado por arbustos de S. aste- 
racae, P. gentianoides , Lupinus montanus y árboles de Alnusjorullensis, P. hartwegii y Abies religiosa. 
Los histogramas de alturas definieron dos y tres estratos verticales, respectivamente. 

PALABRAS CLAVE: 

Abies religiosa, Alnus joruiiensis, estructura vegetal, Pinus hartwegii, Estado de México. 


ABSTRACT 

In the period from February through June 2004, the presence of different bird species was 
studied in two zones located at the Experimental Forest Station Zoquiapan (EFEZ for its Spanish 
acronym), State of México, one with slightly disturbed pine forest (ZOQ 1) and another with mixed 
perturbed forest (ZOQ 2). The points recounting method with fixed radio for bird monitormg was used 
as well as the square embedded method and the quadrant with a central point, the last one to allow 
sampling the vegetation. The distribution of birds in vertical vegetation strata was characterized and 
compared using Poisson regression analysis and frequency observation: furthermore, a specific 
analysis was developed to describe the vegetal structure in both zones. The vertical distribution of 
birds followed four patterns: low (B), high (A), high-low (AB) and ¡ndeterminate (I). The individual’s 


1 Colegio de Postgraduados, km. 36.5, carretera México-Texcoco, Montecillos, 56230, Estado de México. 
Postgrado en Ganadería, Ce para correspondencia: saulu@colpos.mx 

2 Postgrado Forestal. 

3 Postgrado en Estadística. 
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Distribución vertical de aves en un bosque templado ... 


abundance per stratum had a different 
valué within strata and zones; in both 
sites the upper stratum showed a greater 
abundance than in the low and médium 
ones, although the number of individuáis 
per specie registered in each stratum 
was the greatest in the perturbed mixed 
forest. The relative importance valué 
(VIR for its Spanish acronym) of vegeta- 
ron confirmed than ZOQ 1 had Pinus 
hartwegii as the predominant specie with 
some shrubs of Penstemon gentianoides 
and Senecio asteracae. ZOQ 2 was a 
mixed forest comprised by S. asteracae, 

P. gentianoides, Lupinus montanus 
shrubs and Alnus jorullensis, P hartwegii 
and Abies religiosa trees. The height 
histograms determined two and three 
vertical strata, respectively. 

KEY WORDS: 

Abies religiosa, Alnus jorullensis, vege- 
taron structure, Pinus hartwegii, México 
State. 

INTRODUCCIÓN 

México se caracteriza por tener una 
extraordinaria riqueza biótica asociada a 
su cobertura forestal que abarcaba apro- 
ximadamente 120 millones de ha, de las 
cuales 55,3 millones correspondían a 
bosques templados. Actualmente sólo 
quedan alrededor de 21 millones de ha de 
estos bosques con vegetación primaria, 
ya que el resto corresponde a zonas con 
vegetación perturbada (SARH, 1994). 
Entre los principales agentes de perturba- 
ción de los ecosistemas en México se 
encuentran las actividades agropecua- 
rias, la tala inmoderada y los incendios 
naturales e inducidos. Éstos son eventos 
que han provocado un alto decremento 
de la cobertura forestal (Masera, 2001), 
principalmente la de bosques templados 
(INEGI-UNAM, 2001). 

Los remanentes de estos bosques se 
localizan en regiones templadas, subhú- 
medas y frías de las montañas: 1 500 a 
4 000 msnm (Miranda y Hernández-X., 
1963; Rzedowski, 1978); por lo que su 
distribución coincide fisiográficamente con 
las principales cordilleras que atraviesan al 


país (González, 2003). Así, una superficie 
relevante de vegetación templada se loca- 
liza en el Eje Neovolcánico Transversal 
(ENT) (Ferrusquía, 1998; González, 2003); 
región que alberga al sistema montañoso 
conocido como Sierra Nevada, sección que 
ha sido altamente perturbada por fenó- 
menos naturales y actividades antrópicas 
como las ya antes señaladas. Por tanto, sus 
bosques constituyen sistemas idóneos para 
el estudio de perturbaciones y sus efectos 
sobre las aves, ya que presentan áreas 
usualmente en forma de parches de vege- 
tación cuya composición florística permite 
identificar diversos estados sucesionales. 
Tal es el caso de los bosques mixtos de aile 
( Alnus jorullensis), pino ( Pinus hartwegii) y 
oyamel {Abies religiosa) que representan 
un estado sucesional intermedio, el cual 
contrasta con el bosque de pino (P. hart- 
wegii) con zacatonal {Muhlenbergia 
macroura), que representa un estado suce- 
sional tardío de mayor extensión en dicha 
sierra (Zavala, 1984; Rzedowski, 1996). Sin 
embargo, no se han realizado estudios en la 
zona de interés referentes a la forma en que 
las comunidades de aves responden a 
modificaciones producidas por perturba- 
ciones en la estructura y composición de la 
vegetación. 

Los disturbios, según Pickett y White 
(1985), son procesos ecológicos dispara- 
dores de la sucesión secundaria, los 
cuales producen mosaicos de vegetación 
con diferente complejidad estructural (i.e., 
estados sucesionales); a su vez, estos 
parches juegan un papel importante en el 
mantenimiento de la biodiversidad, tanto 
a nivel local como regional. Las modifica- 
ciones en la estructura de la vegetación 
provocan cambios en las poblaciones de 
algunas especies animales, las cuales 
pueden verse beneficiadas, afectadas 
negativamente e incluso no exhibir alguna 
respuesta aparente (Marsden, 1998; 
Flores et al., 2002). La manera específica 
en que las perturbaciones naturales o 
inducidas afectan la ecología de las aves 
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ha sido en general poco investigada, 
menos aún sus implicaciones para la 
conservación de este grupo taxonómico. 

Por lo anterior, es imperativo que los 
estudios de las relaciones ave-hábitat 
deban incluir análisis de impacto- 
respuesta, así como de los mecanismos 
subyacentes a las modificaciones de la 
estructura del hábitat producidas por 
perturbaciones (Wiens, 1989). Esto se 
fundamenta en la estrecha relación que 
existe entre las aves, la estructura de la 
vegetación y los recursos asociados 
disponibles (por ejemplo alimento) antes 
y después de las perturbaciones (Brush y 
Stiles, 1986; Holmes, 1990). Por ello, 
estudios de este tipo deben considerar la 
composición y estructura vegetal 
(Shugart y Patten, 1972), las cuales están 
íntimamente ligadas con la disponibilidad 
de recursos alimenticios y con la hetero- 
geneidad ambiental que en conjunto 
determinan la distribución de las aves en 
un plano horizontal, vertical o temporal 
(Cody, 1968; Karr, 1971). Por lo tanto, es 
primordial estudiar el uso de hábitat y la 
distribución vertical de este grupo taxonó- 
mico en condiciones de perturbación 
(Collar et al., 1994; Stotz etal., 1996;Alta- 
mirano et al., 2002). 

OBJETIVO 

Caracterizar y contrastar la distribu- 
ción vertical de aves, respecto a la propor- 
ción de individuos observados en los 
perfiles de estratificación de la vegetación, 
bajo dos condiciones: bosque mixto (pertur- 
bado) y bosque de pino (poco alterado). 

METODOLOGÍA 

El área de estudio se ubica en la 
porción sur del Parque Nacional 
Zoquiapan y Anexas, entre las coorde- 
nadas geográficas de 19°12'30” a 
19°20'00” de latitud norte y de 98°30'00” 


a 98°42'30” de longitud oeste, que 
corresponden a los límites del Estado de 
México y Puebla, a una altitud que va de 
3 080 a 3 690 msnm. Forma parte del 
ENT, específicamente del macizo monta- 
ñoso denominado Sierra Nevada (Blanco 
etal., 1981). En dicha área, dos zonas de 
22 ha cada una fueron delimitadas con 
base en la composición de su cubierta 
arbórea y perturbación aparente: a) 
bosque de pino [Pinus hartwegii (Lindl.)] 
con poca alteración en su estructura 
(ZOQ 1) y b) bosque mixto [Alnus joru- 
II ensis (H.B.K.) - P. hartwegii - Abies reli- 
giosa (H.B.K.; Cham y Schl.)] con vegeta- 
ción secundaria producto de perturbación 
tanto natural como humana (ZOQ 2). 

Comunidades vegetales 

Para el muestreo de vegetación 
herbácea y arbustiva, en cada zona se 
utilizó la metodología de cuadros empo- 
trados (Oosting, 1956) y para arbóreas, 
cuadrantes con punto central (Cottam y 
Curtís, 1956), por ser éstas las más 
adecuadas para evaluaciones forestales 
en bosques templados (Zavala, 1984). 
Mediante recorridos preliminares en 
campo se definió un universo de 21 
unidades de muestreo (UM) circulares por 
zona distribuidas sistemáticamente, las 
cuales tuvieron un radio de 37,5 m y una 
superficie de 0,45 ha, donde se regis- 
traron e identificaron taxonómicamente 
todos los árboles, arbustos y herbáceas. 
La toma de muestras vegetales se realizó 
una sola vez durante el mes de mayo del 
2004, identificando las especies en los 
herbarios de la Universidad Autónoma 
Chapingo (CHAP) y del Colegio de 
Postgraduados (CHAPA). 

Comunidades de aves 

Los datos de distribución y abun- 
dancia de aves dentro de las UM emplea- 
das para la vegetación, se obtuvieron 
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Distribución vertical de aves en un bosque templado ... 


mediante la metodología de recuento en 
puntos con radio fijo de 25 m (Hutto et al., 
1986), por ser ésta la de mayor espectro 
de muestreo en la observación de aves 
en un plano vertical y horizontal, sobre 
todo a nivel del dosel (Rappole et al., 
1998). Se establecieron un total de 12 
puntos de recuento (PR) por zona, distri- 
buidos sistemáticamente a manera de 
cuadrícula a distancias de 200 m entre 
cada PR, con la finalidad de obtener una 
cobertura equitativa en toda el área y una 
adecuada independencia estadística en 
los datos (Ralph et al., 1996). Del 14 de 
febrero al 19 de junio del 2004 se reali- 
zaron un total de 12 muéstreos por zona, 
durante los cuales se cambió el orden en 
que los puntos fueron visitados; es decir, 
durante el primer muestreo se inició con 
el punto uno, en el segundo muestreo con 
el punto dos, y así sucesivamente, con el 
fin de disminuir los efectos del sesgo por 
secuencia-horario (Gram y Faaborg, 
1997). Las observaciones fueron reali- 
zadas con binoculares Bushnell 25 x 50 y 
la identificación de especies con guías de 
campo estándar (i.e., Peterson y Chalif, 
1989; Howell y Webb, 1995; National 
Geographic Society, 2002). 


ANÁLISIS DE DATOS 
Comunidades vegetales 

Para caracterizar la fisonomía 
vegetal, en cada zona se consideraron 
para árboles y arbustos, los criterios utili- 
zados en el Centro Agronómico Tropical 
de Investigación y Enseñanza; el cual 
considera a las siguientes variables y 
análisis como los más importantes para 
dicho fin en bosques templados (Beer, 
1984): diámetro medido a la altura del 
pecho (DAP=1,3 m sobre el nivel del 
suelo), área basal calculada a partir de la 
variable anterior 

ab = — DAP 2 
4 

y altura total medida utilizando un altí- 
metro tipo Haga (Greig-Smith, 1983). 

Los valores de densidad, domi- 
nancia y frecuencia, así como de impor- 
tancia por especie se determinaron como 
abajo se presentan (Stiling, 1999; Villavi- 
cencio-Enríquez y Valdez-Hernández, 
2003): 


s 

Densidad = Número de individuos / Area muestreada 

Densidad relativa = (Densidad por especie / Densidad de todas las especies) * 100 

✓ 

Dominancia = Total el área basal / Area muestreada 
Dominancia relativa = (Dominancia por especie / Dominancia de todas las especies) * 100 
Frecuencia = UM en que está presente la especie / Número total de UM 
Frecuencia relativa = (Frecuencia por especie / Frecuencia de todas las especies) * 100 


Valor de importancia = Densidad relativa + Frecuencia relativa 
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La ecuación del valor de importancia 
proporciona valores con relación a 300 %, 
para manejar aquéllos con referencia a 
100 % se utilizó la siguiente ecuación: 

Valor de importancia = 

(Densidad relativa + Dominancia relativa 
+ Frecuencia relativa) / 3 


Asimismo, y con objeto de describir la 
estructura vertical de la vegetación en 
ambas zonas, se realizaron histogramas 
con la altura de los individuos censados 
(Zarco, 2007) y se definieron tres estratos 
(Figura 1 ): a) inferior (< 3 m), b) medio (de 
3 a 9 m) y c) superior (> 9 m). 



Figura 1. Estratificación vertical de la vegetación en bosque poco alterado (ZOQ 1) y 
con perturbación (ZOQ 2) de la Estación Forestal Experimental Zoquiapan (EFEZ). 
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Distribución vertical de aves en un bosque templado ... 


Comunidades de aves 

La distribución vertical de las aves 
se caracterizó considerando los estratos 
mencionados anteriormente (MacArthur 
y MacArthur, 1961; Rappole et al., 
1998). El número promedio y la propor- 
ción de individuos por especie se esti- 
maron y contrastaron por estrato en 
cada zona mediante modelos de regre- 
sión Poisson, utilizando solamente los 
datos de las especies de aves más 
abundantes ( i.e ., aquéllas con más de 
21 individuos). Posteriormente, con 
estos datos se modeló la abundancia de 
las aves, y ante los estratos verticales 
de la vegetación (x) en cada zona, reali- 
zando un ajuste de la variable “y” a las 
variables “x”, mediante un modelo lineal 
generalizado (GLM). En dicho análisis 
se consideraron como coeficientes esta- 
dísticamente significativos aquellos en 
los que p<0,05 (McCullagh y Nelder, 
1989; González-Oreja, 2003). En todos 
los casos se supuso una distribución 
Poisson para los valores de abundancia 
de especies, con el logaritmo como 
función de relación (vínculo). Además, 
no se predeterminó en ningún caso, la 
estructura del modelo, el cual fue ajus- 
tado mediante un proceso de regresión 
Poisson cuadrática por pasos hacia 
atrás con los tres estratos de la vegeta- 
ción (Guisan y Zimmermann, 2000), 
para determinar cuál de los estratos 
explica mejor la variabilidad en la rela- 
ción: abundancia de aves-estratos de la 
vegetación (Rotenberry y Wiens, 1980). 
Se utilizaron modelos Poisson porque 
tienen la ventaja de analizar datos 
procedentes de experimentos Bernoulli 
con probabilidades de éxito muy 
pequeñas (como los obtenidos en el 
presente estudio); por lo que se les 
conoce como modelos de eventos raros 
(Said y Zárate, 2003). Dicho análisis se 
llevó a cabo utilizando el procedimiento 
GLM del paquete estadístico R-versión 
2.0.1. (2004). 


La ocurrencia de las aves en los 
estratos de la vegetación de cada zona se 
determinó mediante la Frecuencia 
de Observación (FO) de las especies de 
aves abundantes (< 21 individuos) y poco 
abundantes (< 21 individuos); la ecuación 
para realizar los cálculos fue la siguiente 
(Curts, 1993): 

FO = (No. de individuos registrados en un 
estrato / No. total de individuos) * 100 

RESULTADOS 

La composición florística para ZOQ 
1 y ZOQ 2 se muestra en la Tabla 1 . Las 
especies vegetales que tuvieron los 
mayores valores de importancia en la 
zona ZOQ 1 fueron el árbol Pinus hart- 
wegii, seguido por los arbustos Pens- 
temon gentianoides (H.B.K.) ( Poiret.) y 
Senecio asteracae (H.B.K.) (Poiret.) 
(Tabla 2). En contraste, las especies que 
tuvieron los mayores valores de impor- 
tancia en la zona ZOQ 2 fueron los 
arbustos S. asteracae, P gentianoides y 
Lupinus montanus (H. B. K.), seguidos 
por los árboles A. jorullensis, P hartwegii 
y A. religiosa (Tabla 3). En la Tabla 4 se 
indican los componentes del modelo: 
estrato, zona, interacción (estrato*zona) 
que tuvieron efecto sobre la distribución 
vertical de las aves. 

Los resultados permiten señalar que 
la distribución de las aves en los estratos 
verticales de la vegetación siguió cuatro 
patrones de comportamiento en las dos 
zonas de estudio. En el primero (patrón 
bajo “B”), sólo tres especies de aves se 
encuentran en mayor proporción en el 
estrato inferior, seguido por el medio y, 
por último, el superior (Figura 2). En el 
segundo (patrón alto “A"), ocho especies 
de aves se encuentran en mayor propor- 
ción en el estrato superior, seguido por el 
medio y, por último, el inferior (Figura 3). 
En el tercero (patrón alto y bajo “AB”), 
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Tabla 1. Composición florística de estratos verticales en bosque poco alterado (ZOQ 1) 
y con perturbación (ZOQ 2) de la Estación Forestal Experimental Zoquiapan (EFEZ). 

Familia Especie Forma 

de vida ZOQ 1 ZOQ 2 


Estrato superior 

Pinaceae 

Pinus hartwegii (Lindl.) 

Árbol 

X 

X 

Pinaceae 

Abies religiosa (H. B. K.) (Cham & Schl.) 

Árbol 


X 

Betulaceae 

Alnus jorullensis (H. B. K.) 

Árbol 


X 

Estrato medio 

Melanthiaceae 

Stenanthium frigidum (Shltdl. & Cham) (Kunth) 

Arbusto 

X 

X 

Scrophulariaceae 

Castilleja moranensis (H. B. K.) 

Arbusto 

X 


Scrophulariaceae 

Penstemon gentianoides (H. B. K.) (Poiret.) 

Arbusto 

X 

X 

Ericaceae 

Vaccinium caespitosum (Michx) 

Arbusto 

X 


Fabaceae 

Lupinus montanus (H. B. K.) 

Arbusto 

X 


Asteraceae 

Senecio cinerarioides (H. B. K.) 

Arbusto 


X 

Asteraceae 

Senecio asteracae (H. B. K.) (Poiret.) 

Arbusto 

X 

X 

Asteraceae 

Chaptalia ehrenbergia (Sch. Bip.) (Hemsl.) 

Arbusto 


X 

Asteraceae 

Senecio angulifolius (D. C.) 

Arbusto 


X 

Asteraceae 

No determinada 

Arbusto 

X 


Caprifoliaceae 

Symphoricarpos microphyllus (H. B. K.) 

Arbusto 


X 

Estrato inferior 

Geraniceae 

Geranium potentillifolium (D. C.) 

Herbácea 

X 


Geraniceae 

Geranium cruceroense (Kunth) 

Herbácea 


X 

Bryaceae 

Musgo sp. 

Herbácea 

X 

X 

Rosaceae 

Fragaria mexicana (Schlecht) 

Herbácea 

X 

X 

Oxalidaceae 

Oxalis jacquiniana (H. B. K.) 

Herbácea 

X 


Commelinaceae 

Weldenia candida (Schult F.) forma candida 

Herbácea 

X 


Onagraceae 

Oenothera roseae (L'Hérex Ait.) 

Herbácea 

X 

X 

Asteraceae 

Cirsium nivale (H. B. K.) (Sch. Bip.) 

Herbácea 

X 

X 

Umbelliferae 

Tanschia nudicaulis (Schelecht.) 

Herbácea 

X 


Umbelliferae 

Eryngium carlinae (Delar. F.) 

Herbácea 

X 


Fabaceae 

Trifolium sp. 

Herbácea 

X 

X 

Fabaceae 

Trifolium amabile (H. B. K.) 

Herbácea 

X 

X 

Poaceae 

Calamagrostis tolucensis (Kunth) 

Herbácea 

X 


Poaceae 

Muhlenbergia macroura (H. B. K.) (Kunth) 

Herbácea 

X 

X 

Poaceae 

Festuca tolucensis (H. B. K.) 

Herbácea 


X 

Violaceae 

Viola grahami (Benth.) 

Herbácea 

X 


Caryophyllaceae 

Stellaría crespidata (Willd) 

Herbácea 


X 

Labiatae 

No determinado 

Herbácea 


X 

Solanaceae 

Solanum demissum (Lindl.) 

Herbácea 


X 

Hydrophyllaceae 

Phacelia platycarpa (Cau.)(Spreng.)(E. Guizar) 

Herbácea 


X 

Rubiaceae 

Galium aschenbornii (Schauer.) 

Herbácea 


X 
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Tabla 4. Análisis de efectos (regresión Poisson con modelos GLM) para las 21 especies de aves más 
abundantes en ambas zonas: poco alterada (ZOQ 1) y con perturbación (ZOQ 2) de la Estación Forestal 

Experimental Zoquiapan (EFEZ). 


Especie Desviación de los residuales 

Estrato 

Zona 

Interacción 

Poecite sclateri 

2,34 

0,002785 ** 

- 

- 

Peucedramus taeniatus 

2,74 

0,001210 ** 

1 ,34e-08 *** 

- 

Picoides villosus 

2,22 

2,24e-06 *** 

0,00140 ** 

- 

Junco phaeonotus 

3,86 

2,89e-12 *** 

0,02927 * 

0,00712** 

Ergaticus ruber 

3,35 

0,0114 * 

9,62e-05 *** 

2,64e-11 *** 

Turdus migratoríus 

3,39 

0,036622 * 

0,500777ns 

0,001372 ** 

Certhia americana 

3,46 

0,02389 * 

0,00123 ** 

- 

Oriturus superciliosus 

3,10 

0,00109 ** 

- 

- 

Catharus occidentalis 

3,45 

0,04812 * 

0,00479 ** 

- 

Regulus caléndula 

3,96 

- 

4,43e-05 *** 

- 

Troglodytes aedon 

4,49 

9,20e-14 *** 

80e-05 *** 

- 

Dendroica coro nata 

4,50 

5,68e-05 *** 

1 ,74e-08 *** 

- 

Dendroica townsendi 

4,01 

0,049935 * 

- 

- 

Sialia mexicana 

5,20 

0,0113* 

6,72e-05 *** 

- 

Empidonax spp. 

5,58 

0,00149 ** 

9,14e-09 *** 

- 

Carduelis pinus 

3,66 

0,301ns 

- 

- 

Cyanocitta stelleri 

2,76 

0,995ns 

- 

- 

Myioborus miniatus 

1,78 

0,994ns 

- 

- 

Pipilo erythropthalmus 

3,80 

0,288ns 

- 

- 

Sitta pygmaea 

4,53 

0,00454** 

- 

- 

Vermivora celata 

2,14 

0,994ns 

0,00346 ** 

- 


***p=0,001, **p=0,01, *p=0,05, ns=sin efecto. 


cuatro especies de aves se encuentran 
en mayor proporción en el estrato inferior 
o el superior, pero siempre utilizan el 
estrato medio en menor grado que los 
otros dos (Figura 4). En el cuarto y último 
(patrón indeterminado “I”), seis especies 
de aves se encuentran de manera similar 
en los tres estratos de la vegetación en 
ambas zonas (Figura 5). 

En general, los resultados de regre- 
sión Poisson con modelos GLM no 
evidenciaron alguna tendencia numérica 
en la presencia de aves en los estratos 
entre zonas. Sin embargo, en ambas 
condiciones del bosque, poco alterado y 
perturbado, las aves siguieron los cuatro 
patrones de distribución mencionados 


anteriormente. Además, la proporción de 
individuos en los estratos al interior de 
cada zona fue similar: mayor presencia 
de aves en el estrato superior que en el 
inferior y medio (Tabla 5, Figura 6). La 
Tabla 6 muestra las diferencias significa- 
tivas entre estratos y zonas para cada 
especie de ave registrada como más 
abundante. Los resultados obtenidos 
mediante frecuencia de observación 
(Figura 7) muestran concordancia sola- 
mente con aquellos generados usando 
modelos GLM (Figura 6) para las propor- 
ciones de individuos en el estrato inferior; 
en los estratos medio y superior las 
proporciones obtenidas fueron diferentes 
entre ambos métodos. 






NCincro Estimado de Incivitkns 


4 


Distribución vertical de aves en un bosque templado 



Figura 2. Patrón B de distribución vertical de la abundancia de aves en bosque poco 
alterado (ZOQ 1 —I—) y con perturbación (ZOQ 2 -4-) de la Estación Forestal 
Experimental Zoquiapan (EFEZ): a) Troglodytes aedon, b) Oriturus superciliosus 

y c) Junco phaeonotus. 
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Figura 3. Patrón A de distribución vertical de la abundancia de aves en bosque poco 
alterado (ZOQ 1 -4- ) y con perturbación (ZOQ 2 —I— ) de la Estación Forestal 
Experimental Zoquiapan (EFEZ): a) Poecile sclateri, b) Sialia mexicana , c) Empidonax 
spp., d) Dendroica coronata, e) Peucedramus taeniatus, f) Dendroica townsendi, 

g) Ergaticus ruber y h) Sitia pygmaea. 
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Distribución vertical de aves en un bosque templado 



Figura 4. Patrón AB de distribución vertical de la abundancia de aves en bosque poco 
alterado (ZOQ 1 -4- ) y con perturbación (ZOQ 2 -4— ) de la Estación Forestal 
Experimental Zoquiapan (EFEZ): a) Certhia americana, b) Picoides villosus, 
c) Catharus occidentales y d) Turdus migratorius. 
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Figura 5. Patrón I de distribución vertical de la abundancia de aves en bosque poco 
alterado (ZOQ 1 -I—) y con perturbación (ZOQ 2 -t— ) de la Estación Forestal 
Experimental Zoquiapan (EFEZ): a) Regulus caléndula, b) Pipilo erythropthalmus, 
c) Cyanocitta stelleri, d) Carduelis pinus, e) Myioborus miniatus y f) Vermivora celata. 
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Distribución vertical de aves en un bosque templado ... 


Tabla 5. Proporción de individuos por estrato vertical obtenida con GLM para las 21 
especies de aves más abundantes en bosque poco alterado (ZOQ 1) y con 
perturbación (ZOQ 2), de la Estación Forestal Experimental Zoquiapan (EFEZ). 


ESPECIE 

Inferior 

ZOQ 1 
Medio 

Superior 

Inferior 

ZOQ 2 
Medio 

Superior 

Poecile sclateri 

0,2 

0,9 

5,7 

0,2 

0,9 

5,7 

Peucedramus taeniatus 

0,1 

2,0 

3,5 

0,0 

0,2 

0,4 

Picoides villosus 

0,1 

0,0 

1,1 

0,2 

0,1 

2,9 

Junco phaeonotus 

6,2 

1,7 

0,8 

8,6 

2,0 

0,2 

Ergaticus ruber 

0,1 

0,1 

0,8 

1,5 

5,7 

8,5 

Turdus migratorius 

2,1 

1,0 

2,6 

2,5 

0,2 

0,7 

Certhia americana 

1,2 

0,5 

3,5 

0,6 

0,2 

1,7 

Oriturus superciliosus 

2,6 

1,7 

0,8 

0,0 

0,0 

0,0 

Catharus occidentales 

0,3 

0,1 

0,2 

0,9 

0,3 

0,6 

Regulus caléndula 

1,0 

1,0 

1,0 

0,1 

0,1 

0,1 

Troglodytes aedon 

4,4 

1,2 

0,1 

8,0 

2,3 

0,3 

Dendroica coronata 

0,1 

2,1 

4,0 

0,0 

0,4 

0,8 

Dendroica townsendi 

0,0 

0,3 

1,6 

0,0 

0,3 

1,6 

Sialia mexicana 

0,4 

1,0 

1,3 

0,0 

0,1 

0,2 

Empidonax spp . 

0,0 

0,3 

1,1 

0,1 

1,4 

5,0 

Carduelos pinus 

0 

0 

0 

2,8 

2,8 

2,8 

Cyanocitta stellerl 

1,8 

1,8 

1,8 

1,3 

1,3 

1,3 

Myioborus miniatus 

0 

0 

0 

2,4 

2,4 

2,4 

Pipilo erythropthalmus 

0 

0 

0 

5 

5 

5 

Sitta pygmaea 

0,1 

0,9 

1,5 

0,0 

0,0 

0,0 

Vermivora celata 

0,1 

0,1 

0,1 

1,7 

1,7 

1,7 

TOTAL 

20,8 

16,7 

31,5 

35,9 

27,4 

41,9 

% 

29,9 

24,3 

45,7 

34,2 

26,1 

39,7 
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Figura 6. Distribución vertical de la abundancia de aves en bosque poco alterado (ZOQ 1 ) y 
con perturbación (ZOQ 2) de la Estación Forestal Experimental Zoquiapan (EFEZ). 
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Tabla 6. Análisis de diferencias significativas entre estratos (inferior, medio, superior) y 
zonas [poco alterada (ZOQ 1), perturbada (ZOQ 2)] para las 21 especies de aves más 
abundantes (regresión Poisson con modelos GLM) de la Estación Forestal Experimental 

Zoquiapan (EFEZ). 


Especie 

Intercepto 

inferior vs medio 

Medio vs superior 

Poco alterada 
vs perturbada 

Poecile sclateri 

0,000452*** 

0,002785** 

4,04e-13*** 

. 

Peucedramus taeniatus 

0,009442** 

0,001210** 

0,000139*** 

1 ,34e-08*** 

Picoides villosus 

1,16e-05**‘ 

0,21500ns 

2,24e-06*** 

0,00140** 

Junco phaeonotus 

<2e-16*** 

2,30e-11*** 

2,89e-12*** 

0,02927* 

Ergaticus ruber 

0,3173ns 

1,0000ns 

0,0114* 

9,62e-05*** 

Turdus migratorius 

0,000243*** 

0,036622* 

0,423572ns 

0,500777ns 

Certhia americana 

0,36026ns 

0,02389* 

2,15e-05*** 

0,00123** 

Oríturus superciliosus 

1 ,26e-07*** 

0,12650ns 

0,00109** 

- 

Catharus occidentalis 

0,00206** 

0,04812* 

0,30107ns 

0,00479** 

Regulus caléndula 

0,868ns 

4,43e-05*** 

- 

- 

Troglodytes aedon 

<2e-16*** 

5,82e-13*** 

9,20e-14*** 

5,80e-05*** 

Dendroica coro nata 

0,00767** 

0,00082*** 

5,68e-05*** 

1 ,74e-08*** 

Dendroica townsendi 

0,001483** 

0,049935* 

0,000297*** 

- 

Sialia mexicana 

0,0281* 

0,0481* 

0,0113* 

6,72e-05*** 

Empidonax spp. 

1 ,96e-06*** 

0,00149** 

5,19e-07*** 

9,14e-09*** 

Carduells pinus 

0,388ns 

0,819ns 

0,30 Ins 

- 

Cyanocitta stelleri 

0,0135* 

0,3196ns 

- 

- 

Myioborus minia tus 

0,994ns 

0,995ns 

0,994ns 

- 

Pipilo erythropthalmus 

0,0130* 

0,9951ns 

0,0134* 

0,3414ns 

Sitia pygmaea 

0,01296* 

0,02169* 

0,00454** 

- 

Vermivora celata 

0,994ns 

0,994ns 

0,994ns 

- 


***p=0,001, **p=0,01, *p=0,05, ns=sin efecto. 



Figura 7. Porcentaje de abundancia de aves por estrato mediante frecuencia de 
observación (FO) en bosque poco alterado (ZOQ 1) y con perturbación (ZOQ 2) de la 
Estación Forestal Experimental Zoquiapan (EFEZ). 
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DISCUSIÓN 

Lentijo y Kattan (2005) evaluaron la 
relación entre la estructura de la vegeta- 
ción y la estratificación vertical de la 
avifauna en parches de bosque con rege- 
neración natural y bosque maduro, en la 
cordillera central de Colombia. Encon- 
traron que los bosques difirieron en 
aspectos estructurales, principalmente en 
la cobertura del dosel. Sin embargo, no 
encontraron diferencias en la forma de 
cómo las aves usaban los dos hábitats. 
Sus resultados indican que los bosques 
maduros y con regeneración son funcio- 
nalmente equivalentes en cuanto a la 
estratificación de las aves. Con esta infor- 
mación señalan que los segundos no 
necesariamente tienen un efecto negativo 
sobre las comunidades de aves ni sobre 
sus procesos ecológicos. Sin embargo, 
diferimos de dicha afirmación ya que 
nuestras evidencias y resultados pueden 
ser mejor explicados con base en la 
complejidad de la estructura de la vegeta- 
ción, la cual determina la distribución y 
dinámica de las aves en las dimensiones 
vertical y horizontal (MacArthur, 1965). 

Estos resultados aportan informa- 
ción relevante de cómo en el bosque 
poco alterado (ZOQ 1), ambas dimen- 
siones tienden a ser más homogéneas 
que en el bosque perturbado (ZOQ 2), 
originando que la abundancia de indivi- 
duos en los diferentes estratos verticales 
de la vegetación sea menor en el primero 
con respecto al segundo. Greenberg 
(1981) argumenta que bosques tropicales 
maduros sin alteración exhiben un mayor 
número de individuos de diferentes espe- 
cies de aves en el dosel, mientras que en 
los estratos bajos se presentan menos 
individuos y especies. En contraste, hábi- 
tats con un grado de alteración moderada 
permiten un mayor número de individuos 
de diferentes especies y presentan una 
distribución más homogénea de éstas 
sobre los diferentes estratos; es decir, la 


mayoría de las aves en este tipo de 
bosques generalmente no se limitan a 
algún estrato en particular, sino que se 
vuelven más flexibles en la utilización de 
ellos (García et al., 1998; Bojorges y 
López, 2001). 

En el presente estudio, la zona con 
evidencias de una perturbación pasada 
(posiblemente un incendio) tuvo una 
mayor heterogeneidad en su estructura 
(ZOQ 2); la presencia de vegetación 
secundaria proporcionó un mosaico de 
microhábitats que atrajo algunas espe- 
cies y excluyó a otras pero que, en 
general, mostró un mayor número de indi- 
viduos de aves, quizás debido en parte a 
una mayor cantidad de recursos, particu- 
larmente alimenticios. Así, las perturba- 
ciones modifican la composición y diná- 
mica de una comunidad biótica original 
(Liu et al., 1998), reconociendo una rela- 
ción directa de la biodiversidad con los 
patrones de perturbación recibidos (Lusk, 
1996). Algunos autores (MacArthur y 
MacArthur, 1961; MacArthur et al., 1962) 
han demostrado que la heterogeneidad 
ejerce una importante influencia sobre las 
comunidades de aves, ya que determina 
su distribución tanto horizontal como 
vertical, y han estudiado la selección de 
hábitat utilizando medidas tales como 
composición, densidad y estructura 
vegetal en los planos vertical y horizontal. 

Las aves son generalmente muy 
susceptibles a los efectos de las perturba- 
ciones, tanto a nivel individual (cambios 
en conducta y fisiología) como pobla- 
cional y comunitario (cambios en riqueza, 
diversidad, abundancia y distribución), 
más aún en el caso de especies endé- 
micas o raras, que pueden llegar a 
desaparecer. En consecuencia, se 
producen cambios a diferentes niveles en 
las comunidades de aves, posteriores a la 
perturbación de un bosque (MacFaden y 
Capen, 2002); por ello, las aves pueden 
ser utilizadas como indicadoras del 
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estado que guardan las comunidades 
bióticas (Hill et al., 1997). En general, los 
análisis de las relaciones entre las aves y 
sus hábitats se han centrado en las varia- 
ciones de la respuesta animal en función 
de un conjunto de características del 
hábitat (Willson, 1974). Entre los factores 
que más influyen en los procesos de 
selección de hábitat o nichos específicos 
se encuentran las medidas descriptoras 
de la estructura de la vegetación (Roten- 
berry y Wiens, 1980; Stiles, 1980; Cody, 
1985) y en menor medida la composición 
de especies de plantas (Holmes y 
Robinson, 1981), ya que la distribución y 
disponibilidad de recursos están íntima- 
mente ligados a estas características del 
hábitat (Nocedal, 1984). 

Así pues, la estructura de la vegeta- 
ción es importante en la determinación de 
la estructura y composición de las comu- 
nidades de aves mediante la disponibi- 
lidad de recursos alimenticios (Lambert, 
1992). De esta manera, el impacto de las 
perturbaciones sobre las comunidades de 
aves está asociado al tiempo en que 
ocurren las mismas y a modificaciones en 
la estructura del hábitat, las cuales 
cambian la disponibilidad de recursos. 

Adicionalmente, el número de 
estratos de la vegetación considerado es 
un aspecto determinante en este tipo de 
análisis, debido a que puede alterar los 
resultados. Por ejemplo, Nocedal (1984) 
consideró siete estratos verticales y 
agrupó a las aves en gremios alimenti- 
cios, a diferencia del presente estudio 
que consideró tres estratos verticales y se 
basó en el número de individuos por 
especie, así como en la frecuencia de 
observación de las aves en cada estrato. 
Además, la estructura y composición 
florista del hábitat en la zona perturbada 
(ZOQ 2), condicionó una mayor abun- 
dancia de aves en cada estrato vertical de 
la vegetación, en comparación con la zona 
poco alterada (ZOQ 1) quizás debido a 


una mayor disponibilidad de recursos. 
Como ejemplo, aunque Ergaticus ruber 
se registró en los tres estratos de la vege- 
tación en ambas zonas, en ZOQ 2 se 
encontró en mayor proporción (estratos 
superior y medio), lo cual indica que a 
mayor heterogeneidad vertical del hábitat 
puede existir una mayor variedad de 
recursos alimenticios dispuestos de 
manera más regular y que se encuentran 
disponibles para ser aprovechados por 
las diferentes especies de aves (Pickett y 
Thompson, 1978; Santos y Tellería, 
1997). 

Bajo este contexto, Wunderle y Latta 
(1996) señalan que bosques naturales 
albergan especies de aves especialistas 
de bosques maduros y ciertas especies 
de aves endémicas. En contraste, zonas 
perturbadas proveen de hábitat temporal 
o parcial a muchas especies de aves de 
hábitos flexibles en la utilización de la 
estructura del bosque en regiones sujetas 
a deforestación extrema (Ewel y Putz, 
2004), como es el caso de las zonas bajo 
estudio, particularmente la zona ZOQ 2. 

CONCLUSIÓN 

La distribución vertical de las aves 
en ambas zonas (ZOQ 1 y ZOQ 2) siguió 
cuatro patrones diferentes; pero con 
comportamiento similar entre zonas: bajo 
(B), alto (A), alto-bajo (AB) e indetermi- 
nado (I). El estrato superior tuvo una 
mayor abundancia que el inferior y medio, 
en ambos casos; sin embargo, el número 
de individuos por especie registrado en 
cada estrato fue mayor en el bosque 
perturbado (ZOQ 2). Variaciones en la 
abundancia de individuos por estrato 
entre zonas, responden a cambios en la 
estructura vegetal del hábitat en un plano 
vertical. 

Existe la presencia en el área de 
estudio de dos grupos de aves a nivel de 
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comunidades: a) aquéllas que son muy 
flexibles en su distribución sobre la vege- 
tación en un plano vertical, por lo tanto, se 
localizan en zonas con cierto nivel de 
perturbación, ya que en éstos encuentran 
características parciales o permanentes 
para desarrollarse y; b) especies propias 
de hábitats maduros o en estados clímax 
de la vegetación. 

La distribución vertical exhibida, 
sugiere un papel relevante de las áreas 
perturbadas en la dinámica existente 
entre aves y vegetación a lo largo de los 
procesos ecológicos que conllevan a 
éstas hacia el estado clímax durante la 
sucesión vegetal, debido a que sirven 
como hábitat temporal o permanente para 
muchas especies de aves. Algunas de las 
cuales, son endémicas o se encuentran 
en alguna categoría de riesgo de acuerdo 
con la SEMARNAT (2001; Nom-059- 
Ecol.); tal es el caso de Troglodytes 
aedon. 

Con base en los resultados obte- 
nidos, es importante detectar y evaluar 
otras zonas perturbadas en bosques 
templados del área de estudio, en parti- 
cular, y de la Sierra Nevada, en general, 
con la finalidad de aportar conocimiento 
que contribuya a su mejor manejo y 
conservación. 
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Aplicación de técnicas multicriterio en el 
manejo integral forestal en Durango, México 

Application of multi-criterion techniques on integral forest 

management in Durango, México 

Gustavo Pérez-Verdín 1 *, José Ciro Hernández-Díaz 2 , 
Marco Antonio Márquez-Linares 1 y Aregai Tecle 3 

RESUMEN 

Tradicionalmente, la capacidad del bosque para producir bienes y servicios se ha enfocado 
principalmente a maximizar los beneficios económicos únicamente de la producción de madera. Este 
enfoque no toma en cuenta otras funciones del bosque, ni las interacciones del ecosistema, provo- 
cando eventualmente la degradación de algunos de sus componentes. En este estudio se propone 
una metodología para analizar y resolver un problema forestal que involucra el manejo de varios 
bienes y servicios del bosque, tomando en cuenta la opinión de productores y técnicos forestales. Se 
usaron datos de inventario forestal de un ejido del estado de Durango y se recurrió a la técnica de 
consulta a expertos para construir funciones de respuesta para diez objetivos. Se utilizaron tres varia- 
bles independientes: área basal, número de árboles por hectárea y diámetro cuadrático medio, para 
evaluar cada uno de estos objetivos. Se usó la técnica multicriterio de programación por compromisos 
para determinar la mejor combinación de esas variables que respondan al manejo simultáneo de 
dichos objetivos. Mediante un simulador silvícola y sistemas de información geográfica se evaluaron 
las soluciones considerando la dinámica de crecimiento de los bosques. Las estructuras forestales 
encontradas normalmente fueron superiores a las condiciones actuales de los rodales. Estas diferen- 
cias podrían deberse a prácticas de manejo anteriores o a otros disturbios ambientales que han redu- 
cido los inventarios y la calidad de las masas forestales. El uso de este tipo de herramientas puede 
ayudar a entender la dinámica de los bosques y el manejo simultáneo de varios usos y servicios del 
bosque en el largo plazo. 

PALABRAS CLAVE: 

Dinámica forestal, estructura forestal deseada, Los Altares, programación por compromisos, toma de 
decisiones multicriterio. 


ABSTRACT 

Traditionally, the forest capacity to produce múltiple goods and Services has been overlooked in 
favor of developing alternatives that maximize timber production or economic profit alone. This 
approach is suboptimal, it fails to address ecosystem interactions, and occasionally leads to degrada- 
ron of some of its components. In this study, a technical framework to analyze and solve a multiobjec- 
tive forest problem in an ejido of Durango, México is presented. Forest inventory data and expert 
consultation were used to construct ten objective response functions expressed in terms of tree basal 
area, number of trees per hectare, and quadratic mean diameter. Compromise programming, a 
multiobjective decision-making technique was used to determine the most preferred solution of forest 
structures. Forest vegetation simulator and geographical information systems were used in order to 
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evalúate the Solutions for various clus- 
ters. The simulated optimal forest struc- 
tures were typically above current stand 
conditions. These differences can be due 
to past management practices and other 
forest disturbances that reduced stoc- 
king and quality of residual trees. The 
use of this type of multiobjective decisión 
making approach can help understand 
forest stand dynamics and long-term 
management of forest goods and 
Services. 

KEYWORDS: 

Forest stand dynamics, desired forest 
structures. Los Altares, compromise 
programming, multiobjective decision- 
making 

INTRODUCCIÓN 

Tradicionalmente, la capacidad del 
bosque para producir múltiples bienes y 
servicios se ha enfocado más hacia la 
producción de madera (Mendoza- 
Briseño, 1993), que al corto plazo propor- 
ciona altos beneficios económicos, pero a 
largo plazo causa un detrimento a los 
recursos asociados al bosque. Este 
enfoque no considera las interacciones 
presentes en los ecosistemas forestales y 
eventualmente provoca la degradación de 
algunos de sus componentes. Los 
esquemas de manejo en la década de los 
1960’s ejercieron excesiva presión sobre 
los recursos forestales, hasta el punto de 
que muchas áreas se caracterizaron por 
tener estructuras y especies no desea- 
bles y a menudo sobreexplotadas 
(Gingrich, 1993). Ese esquema preva- 
leció hasta finales de los 1970’s, cuando 
se introdujeron nuevos métodos para el 
cultivo del bosque, con el propósito de 
usar con más eficiencia la capacidad 
productiva de las áreas forestales (Caba- 
llero, 2000). 

Los nuevos métodos, por ejemplo, 
abordaron el tratamiento de rodales 
jóvenes y sobrepoblados, aplicando acla- 
reos y garantizando una mejor calidad de 
la regeneración natural, al dejar en pie los 


mejores árboles semilleros (Caballero, 
2000). Sin embargo, la capacidad del 
bosque para proveer bienes y servicios 
múltiples, incluyendo productos no made- 
rables, no se reconoció en México hasta 
hace algunos años. Los cambios 
recientes en la administración federal y 
local, motivados parcialmente por 
esfuerzos internacionales relacionados 
con el concepto de sustentabilidad en el 
uso y manejo de los recursos naturales, 
generaron nuevos programas institucio- 
nales, como el de Servicios Hidrológicos, 
el de Servicios Ambientales y el de Biodi- 
versidad, en un intento de ayudar a los 
dueños y poseedores del recurso a incre- 
mentar el número de fuentes alternativas 
de ingresos y empleos (Muñoz-Piña et al., 
2008). Estos programas incluyen el forta- 
lecimiento de bienes y servicios, como la 
recreación, ecoturismo, manejo de la 
fauna silvestre, captura de carbono 
atmosférico, producción de agua y biodi- 
versidad, así como la capacitación a los 
dueños del recurso forestal y a los presta- 
dores de servicios técnicos y el reforza- 
miento a la investigación y asistencia 
técnica (CONAFOR, 2000; SEMARNAT, 
2003). 

Recientemente se ha mostrado más 
interés en desarrollar planes o líneas de 
acción estratégicas que contribuyan a 
atender la complejidad del manejo de 
ecosistemas forestales, de tal forma que 
se logren varias metas sociales y ambien- 
tales en forma simultánea, apartándose 
del modelo que se basaba casi exclusiva- 
mente en el crecimiento económico y 
reconociendo la necesidad de integrar 
completamente la protección ambiental y 
el manejo sano de los recursos naturales 
en todas las actividades (Kazana et al., 
2005). Este enfoque holístico debe 
mejorar la capacidad productiva de los 
recursos y diversificar la economía 
regional. Uno de los primeros enfoques 
multicriterio en el manejo forestal en 
México fue presentado por Pérez-Verdín 
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y Tecle (2002), quienes determinaron el 
área basal que podría satisfacer el 
manejo simultáneo de cinco objetivos de 
manejo forestal. Este estudio, basado en 
técnicas no paramétricas, consideró una 
sola variable de decisión y no se inclu- 
yeron simulaciones a través del tiempo. 

En el presente trabajo se usó al ejido 
Los Altares como caso de estudio para 
formular y resolver un enfoque de manejo 
forestal con varios objetivos, conside- 
rando la dinámica del recurso. Este ejido 
se localiza en el noroeste del estado de 
Durango y posee alrededor de 17 800 ha 
arboladas, de las cuales 93 por ciento 
está cubierto con bosques templados, 
principalmente de pino-encino. Las espe- 
cies principales son: Pinos arizonica, P. 
durangensis y P. teocote, que se mezclan 
con vegetación caducifolia como Quercus 
spp. y Populus tremuloides, entre otros 
(Figura 1). Los bienes de este ejido se 
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distribuyen de manera común entre 108 
ejidatarios. 

Los bosques del ejido se han mane- 
jado por más de 40 años, principalmente 
con fines maderables. El aprovecha- 
miento forestal se ha realizado bajo 
esquemas de manejo de bosque regular 
e irregular, que incluyen tratamientos de 
selección, aclareos y árboles padres. La 
corta anual permisible ha decrecido en 
forma constante desde los 16 000 m3 en 
el año 1985 hasta los 6 000 m3 actuales. 
Esta disminución se explica en parte por 
el uso de métodos más refinados de 
inventario forestal, la exclusión o restric- 
ción de diversas áreas para proteger 
cuerpos de agua, caminos, suelos 
frágiles, hábitat de fauna, etc., y la dismi- 
nución de la productividad por disturbios 
naturales y antropogénicos (UAF 
Santiago Papasquiaro, 2005). 
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Figura 1. Ubicación del Ejido Los Altares, Durango. 
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Los objetivos planteados en este 
trabajo son: (i) proveer y demostrar un 
enfoque integral para formular y resolver 
problemas de manejo de un ecosistema 
forestal que involucre la opinión de los 
productores y técnicos forestales, (i¡) 
evaluar la dinámica del bosque bajo un 
esquema de análisis multicriterio, y (iii) 
determinar estructuras forestales desea- 
bles que satisfagan varios objetivos, 
simultáneamente. 

METODOLOGÍA 

En el estudio se usaron datos del 
último inventario forestal del ejido, la 
consulta a expertos y una técnica multicri- 
terio denominada programación por 
compromisos. En esta sección se 
describen los métodos utilizados, los 
objetivos de manejo detectados, las varia- 
bles de decisión (atributos dasométricos) 
y el uso de un simulador silvícola para 
evaluar la dinámica del recurso. 

Programación por compromisos 

La programación por compromisos 
es una técnica multicriterio basada en 
minimizar la distancia entre el nivel de 
desempeño de dos o más objetivos y un 
punto ideal (Yu, 1973; Zeleny, 1974). Alas 
soluciones más cercanas al punto ideal 
se les considera no dominadas y juntas 
constituyen el juego de soluciones 
compromiso (Goicochea et al., 1982; 
Rodríguez-Cotilla, 2000). Una solución no 
dominada es aquella que es factible para 
un problema multiobjetivo y que, siendo la 
mejor, no existe otra solución factible que 
pueda mejorar un objetivo sin afectar 
negativamente alguno de los otros (Tecle 
et al., 1998; Duckstein y Tecle, 2006). Las 
variaciones del nivel de prioridad entre 
0% y 100% producen un juego de solu- 
ciones compromiso, también conocido 
como ‘frontera eficiente’ o ‘frontera de 
Pareto’ (Yu, 1973). 


La figura 2 muestra la relación de un 
problema multiobjetivo (en este caso, de 
dos objetivos) desde la perspectiva de 
programación por compromisos. Si en el 
proceso de optimización, el objetivo y-¡ 
tiene 100% de prioridad (o peso relativo) 
sobre el objetivo y 2 , la solución al 
problema será el eje horizontal. El obje- 
tivo yi tendrá el mejor desempeño mien- 
tras que el objetivo y 2 tendrá el peor 
desempeño. De manera similar, si el obje- 
tivo y 2 tiene 1 00% más prioridad que el 
objetivo yi, la solución será entonces 
el eje vertical. El punto ideal representa el 
mejor desempeño de los dos objetivos 
involucrados, pero por lo general es impo- 
sible de lograr (Goicochea et al., 1982). 
Por lo regular, un objetivo tiende a ceder 
hacia el otro sin que esto provoque solu- 
ciones dominadas. El problema entonces 
es determinar qué punto, a lo largo de la 
frontera eficiente, resulta ser el más efec- 
tivo para ambos objetivos sin tener una 
utilidad negativa (Yu, 1973). Normal- 
mente, ese punto es el que se encuentra 
más cerca al punto ideal (Zeleny, 1982). 

En problemas multicriterio que invo- 
lucran tres o más objetivos la solución se 
torna más complicada, debido a las difi- 
cultades de expresar gráficamente esa 
relación y a la dificultad para encontrar el 
punto ideal. Además, todos los objetivos 
de manejo tienen tendencias diferentes, 
aunque su dirección sea similar. Por 
ejemplo, la biodiversidad y secuestro de 
carbono tienen una relación positiva y 
directa con el diámetro de los árboles. 
Pero el objetivo secuestro de carbono por 
lo general requiere árboles aún más 
grandes (mayor diámetro) que la biodiver- 
sidad. Otras diferencias son más 
evidentes: la recreación y belleza escé- 
nica se perciben mejor en rodales 
abiertos (Perez-Verdín et al., 2008), mien- 
tras que el control de la erosión y la 
productividad del bosque tienen mejores 
calificaciones en bosques cerrados. Para 
solventar este tipo de conflictos, el uso de 
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Figura 2. Representación del punto ideal y de las soluciones no dominadas para dos 
objetivos (yi , yi)- Los puntos A, B y C son las soluciones no dominadas y juntas 
forman la línea conocida como frontera eficiente. El punto B es el más cercano al punto 
ideal y por lo tanto el que más satisface el conjunto de soluciones óptimas. 


técnicas multicriterio como la programa- 
ción por compromisos puede ayudar 
significativamente. En estos casos, la 
distancia más cercana al punto ideal 
puede calcularse por una familia de 
distancias d p expresada como sigue 
(Zeleny, 1982): 



/ 

p ■ 

d p = min 

*T 

i 

í-y.) 


V, , éf >0, 1<p<oo , 


Y V é ' =1 [1] 


donde p es un parámetro de distancia, y¡ 
es el mejor valor o máximo encontrado 
para el objetivo /, mientras y¡ es el valor 
observado del objetivo i, e I es el número 


de objetivos considerados. El parámetro 
es el nivel de prioridad o importancia 
del objetivo i con relación a los otros. El 
valor p es un indicador de escala y repre- 
senta el nivel de importancia que se le 
otorga a las desviaciones individuales 
entre el punto ideal y la frontera eficiente, 
y depende del analista. El concepto de 
distancia para cualquier valor de p, se 
utiliza aquí como una medida para 
expresar preferencias individuales y no 
como un concepto puro de distancia 
geométrica. La distancia, en este caso, 
representa comparaciones, similitudes o 
proximidad con respecto a coordenadas 
individuales, dimensiones y atributos 
(Tecle, 1992; Srinivasa Raju and Pillai, 
1999). El analista puede considerar 
desde minimizar la suma de esa desvia- 
ción (i.e., tener una compensación 
perfecta entre los objetivos) hasta mini- 
mizar la desviación máxima (i.e., sin 
importar la compensación entre los obje- 
tivos) (Yu, 1973; Tecle et al., 1988). Sin 


32 


Aplicación de técnicas multicriterio en et manejo integral forestal en Durango 


embargo, los valores más comúnmente 
utilizados en la programación por compro- 
misos son cuando p= 1, 2 e a>. Cuando 
p-1, es decir, la mayor distancia finita 
entre dos puntos en sentido geométrico 
(Zeleny, 1982) o comúnmente conocida 
como distancia Manhattan o rectilínea, d p 
resulta ser: 



La mínima distancia entre dos 
puntos es una línea recta; esto se logra 
cuando p= 2 y se representa básicamente 
por la ecuación [1]. Cuando p tiende a o°, 
la ecuación [1] se convierte en un 
problema míni-máx (minimizar la máxima 
distancia) y la solución se determina de la 
manera siguiente: 


dp^ = max {éi(y* -y¡)} V, [3] 


A fin de estandarizar valores, la 
ecuación [1] se puede expresar como 
sigue: 


d =min 

p 



y, - y / ) 


[4] 


donde y*¡* representa el peor o mínimo 
valor (mínimo local) del objetivo /. El 
análisis de sensibilidad de esta técnica se 
basa precisamente en evaluar la desvia- 
ción al punto ideal considerando los 
valores más comunes para p y, además, 
evaluando el peso relativo del objetivo i 
en comparación con el resto de los obje- 
tivos (parámetro X¡) (Goicochea et al., 


1982; Zeleny, 1982; Krcmar et al., 2005). 
Todos los cálculos, incluyendo el análisis 
de sensibilidad, fueron realizados en 
hojas de cálculo de Excel®. 


Objetivos de manejo 

Un objetivo de manejo expresa la 
dirección deseada del cambio con 
respecto a una situación actual (Ducks- 
tein y Tecle, 2006). Por lo general indica 
una de tres posibles orientaciones: maxi- 
mizar, minimizar, o mantener su status 
quo (Tecle, 1992; Duckstein y Tecle, 
2006). En el presente estudio se eligieron 
objetivos de manejo para atender varios 
usos del bosque que se espera se 
produzcan en el sistema forestal ejidal. 
Los objetivos seleccionados, sin orden 
jerárquico, son: (1) Mantener la biodiver- 
sidad o riqueza biológica (biodiv), 

(2) Incrementar el secuestro de carbono y 
provisión de oxígeno (carbón), 

(3) Reducir el riesgo de incendios 
(fuego), (4) Incrementar la producción de 
no maderables (no-mad), (5) Aumentar la 
productividad forestal (prod), (6) Incre- 
mentar la recreación (reo), (7) Incre- 
mentar belleza escénica (belleza), 

(8) Disminuir erosión del suelo (suelo), 

(9) Maximizar beneficios económicos de 
la madera (madera) y (10) Elevar la 
calidad y cantidad del agua (agua). 

Estos objetivos fueron considerados 
tomando en cuenta la Ley Forestal y otros 
documentos normativos (SEMARNAT, 
2003), la opinión de los dueños y otros 
interesados, así como trabajos previos 
relativos al manejo forestal con objetivos 
múltiples (Pérez-Verdín y Tecle, 2002). 
Para evaluar su nivel de importancia, 
cuatro grupos de personas conocedoras 
del ámbito forestal participaron mediante 
encuestas asignando un nivel de prio- 
ridad o peso relativo a cada uno de los 
objetivos (Tabla 1). Los grupos de 
expertos fueron: dueños y/o propietarios 
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Tabla 1. Peso relativo de los objetivos de manejo 1 en el ejido Los Altares, Durango. 


Objetivo 

Especialistas 

Dueños y 
propietarios 

Gobierno y 
ONG 

Todos 

Xi 

Orden 

Xi 

Orden 

h 

Orden 

h 

Orden 

Riqueza en 
biodiversidad 

0,104 

6 

0,104 

4 

0,107 

1 

0,104 

4 

Secuestro de 
carbono 

0,094 

7 

0,102 

6 

0,099 

6 

0,100 

7 

Peligro de 
incendios 

0,108 

1 

0,107 

1 

0,106 

2 

0,108 

1 

Productos no 
maderables 

0,084 

10 

0,081 

10 

0,094 

10 

0,082 

10 

Productividad 
del bosque 

0,107 

3 

0,106 

3 

0,100 

5 

0,106 

3 

Recreación 

0,092 

8 

0,093 

9 

0,095 

9 

0,093 

9 

Belleza 
escén ica 

0,092 

9 

0,095 

8 

0,096 

8 

0,094 

8 

Retención del 
suelo 

0,107 

2 

0,104 

5 

0,103 

3 

0,104 

5 

Venta de la 
madera 

0,107 

4 

0,101 

7 

0,098 

7 

0,102 

6 

Abastecimiento 
de agua 

0,105 

5 

0,107 

2 

0,102 

4 

0,107 

2 


1 La estandarización se hizo de manera que la suma de los pesos relativos sea igual a uno. 


de bosque, especialistas forestales, 
representantes de gobierno y organiza- 
ciones no gubernamentales. En total, se 
entrevistaron 144 personas. 

Consulta a expertos 

Muchos autores sugieren que en 
ausencia de información empírica que 
permita construir expresiones matemá- 
ticas para explicar el comportamiento de 
bienes y servicios del bosque con 
respecto a ciertos atributos dasométricos, 
puede usarse la consulta a expertos (ver 
Alho y Kangas, 1997; Leskinen y Kangas, 
1998; Kangas et al., 2000; Kazana et al., 
2005; Ray Burgman, 2006; Zadnik, 2006). 

La opinión de expertos permite 
generar o respaldar opciones viables en 


una cierta área de interés. La consulta a 
expertos también se ha utilizado para 
determinar los rangos de variables 
ambientales que son apropiados para un 
determinado uso del suelo (Bojórquez- 
Tapia et al., 2001), ordenamiento ecoló- 
gico (Márquez-Linares et al., 2006), o 
para tomar decisiones consensuadas 
cuando existen conflictos por el uso del 
territorio (Bojórquez-Tapia et al., 1994; 
Mendoza et al., 2005). Sin embargo, es 
frecuente preguntarse cómo identificar los 
expertos a ser consultados. En este 
sentido, Kazana et al. (2005, p. 20) 
señalan que los expertos pueden selec- 
cionarse de acuerdo a grado de escola- 
ridad, experiencia, membresías profesio- 
nales, reconocimiento de sus colegas o 
inclusive por autoproclamación. 
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Dada la falta de información empí- 
rica sobre la relación de los bienes y 
servicios ambientales y los atributos 
dasométricos en Durango, se realizaron 
29 entrevistas personales con especia- 
listas forestales que trabajan en la 
entidad. El nivel de educación de los 
especialistas fue desde licenciatura hasta 
doctorado y experiencia profesional 
desde 14 hasta 45 años en el campo del 
manejo forestal. La consulta a expertos, 
junto con la información dasométrica, fue 
usada para determinar funciones de 
producción que eventualmente expre- 
saran estructuras forestales deseables 
para cada uno de los objetivos de 
manejo. 

La pregunta base que se hizo a los 
expertos consultados fue acerca de su 
opinión respecto a cierta condición del 
bosque para maximizar, minimizar o 
mantener cierto objetivo forestal. Para 


ello, se tomaron fotografías de rodales 
con varias condiciones forestales y se 
mostraron las fotos a los expertos para 
obtener sus preferencias y construir una 
ecuación (función de producción) en cada 
objetivo. Las fotografías fueron tomadas 
dentro del ejido tratando de representar la 
diversidad de condiciones forestales 
características del área. Así por ejemplo, 
el área basal varió desde 18,6 m 2 /ha 
hasta 43,9 m 2 /ha; el número de árboles 
fue desde 160 hasta los 2 400 árboles por 
hectárea y el diámetro cuadrático medio 
cubrió una diversidad desde 11.1 cm 
hasta 43.3 cm. Las fotografías editadas 
contenían información dasométrica como 
área basal, número de árboles por 
hectárea y diámetro cuadrático medio 
(Figura 3). Las preferencias de los 
expertos consultados se basaron en una 
escala de 10 puntos donde 1 significa 
muy bajo o muy malo y 10 significa muy 
alto o excelente (Figura 4). 



Figura 3. Ejemplo de una condición base para la construcción de las ecuaciones de 
respuesta. La foto muestra la condición 4 (C4) que cuenta con una área basal de 39,8 
m 2 /ha; número de árboles por hectárea de 270 y un diámetro cuadrático medio de 43,3 cm. 
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Figura 4. Valores de los especialistas forestales para cada uno de los objetivos de 
manejo (A= máximo, B=75% cuartil, C=mediana, D=25% cuartil y E= mínimo. 

El símbolo 0 representa el promedio) 


Variables de decisión 

Las variables de decisión son aque- 
llos atributos dasométricos que pueden 
ser modificadas por el técnico forestal 
para lograr los objetivos de manejo de- 
seados, pudiendo ser discretas o conti- 
nuas (Duckstein y Tecle, 2006). En este 
estudio se consideraron tres atributos de 
rodal: área basal (ba), número de árboles 
por hectárea (nth) y diámetro cuadrático 
medio ( qmd ). Estas tres variables reflejan 
la diversidad de las condiciones ecoló- 
gicas (productividad, densidad y volumen 
del arbolado), son comunes a todos los 
objetivos de manejo y fáciles de manejar 
por el técnico forestal en el campo. Intro- 
ducir más variables al sistema hubiera 
resultado en cierta redundancia dado que 
muchas otras variables dasométricas 
resultan de la interacción de una o más 
de las seleccionadas. Las variables selec- 
cionadas se expresan en m 2 /ha, número 


de árboles por hectárea y centímetros, 
respectivamente. 

Funciones de producción 

A partir de la información recabada 
con la consulta a expertos, se derivaron 
funciones de producción para cada obje- 
tivo del manejo /. Como se mencionó, la 
función de producción y¡ se expresó en 
función de tres variables de decisión 
como sigue: 

y¡ = f¡ (ba, nth, qmd) [5] 

donde, y¡ es la calificación de los expertos 
para el objetivo i, ba es el área basal 
(m 2 /ha), nth es el número de árboles por 
hectárea, y qmd es el diámetro cuadrático 
medio (cm). Se usó la primera derivada 
parcial de cada función y, para determinar 
los valores locales máximos y mínimos, 
usando las expresiones siguientes: 
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y](nth) y] (qmd) - -^L-V , 


dba , 


dnth, 


dqmdi 


[ 6 ] 


dy¡ 


dy¡ 


y¡ (ba) = y, (fiíAi) - y ; (q/77c/) = 


dba 


dnth. 


dy, 

dqmdi 


V ; 


[7] 


La búsqueda de los valores 
máximos o mínimos consiste en encon- 
trar los puntos en los que se presenta un 
cambio en la pendiente a lo largo de la 
función de producción de cada objetivo. 
Dado el carácter multivariado de la expre- 
sión, se requirió un análisis individual de 
cada variable para obtener los valores de 
la ecuación (5). El paquete Solver de 
Excel® fue usado para determinar dichos 
valores para todos los objetivos de 
manejo, sin importar el carácter lineal o 
no lineal de las ecuaciones (Baker, 2006). 
Una vez estimados estos datos, los 
valores resultantes fueron substituidos en 
la ecuación [4] para determinar la 
distancia mínima y eventualmente las 
estructuras forestales deseadas. 


Clasificación y simulación silvícola 

En virtud del alto número de rodales 
que existen en el ejido y con fines prác- 
ticos de manejabilidad, se usó la técnica 
de análisis por conglomerados para 
agrupar rodales con estructuras similares 
en cuanto a número de árboles por 
hectárea, diámetro cuadrático medio, 
área basal, edad y composición de espe- 
cies. Los rodales se agruparon en 15 
conglomerados usando el método de 
relación de promedios por grupos, 
medido por distancias euclidianas 
(Romesburg,1990). En cada conglome- 
rado se identificó un rodal tipo que sirvió 


como representante de las características 
promedio del conglomerado. Cada rodal 
tipo se usó posteriormente en el análisis 
temporal de los datos. La información de 
los rodales tipo incluyó datos generales 
como índices de sitio, índices de 
densidad, pendiente y elevación. Los 
datos a nivel árbol incluyeron especie, 
diámetro y altura de los árboles. 

Se usó el Simulador de Vegetación 
Forestal (FVS) del Servicio Forestal de 
Estados Unidos para simular el creci- 
miento del bosque, la regeneración y la 
mortalidad en los rodales en cada conglo- 
merado. Estas simulaciones permitieron 
conocer las estructuras forestales desea- 
bles a lo largo del tiempo con base en la 
técnica de multicriterio. El FVS consiste 
en una familia de modelos de simulación 
de árboles individuales, independientes 
de la distancia, y puede ser aplicado a 
una amplia variedad de tipos de bosques, 
estructuras de rodales y composición de 
especies (Crookston y Dixon, 2005). Se 
usó la variante de las Montañas Roca- 
llosas del FVS debido a las similitudes de 
los tipos de bosques y a su flexibilidad 
para ajustarse a condiciones específicas. 
Aunque no existen datos comparativos, 
los ecosistemas en esa región y los del 
norte de México comparten caracterís- 
ticas similares en especies, géneros y 
patrones de crecimiento (Richardson, 
1998). Por ejemplo, Pinus arizonica, la 
especie dominante en el norte de 
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Durango, se ha considerado como una 
variedad del P. ponderosa , la especie 
dominante en el suroeste de E.U. (Farjon 
y Styles, 1997; Perry, 1991). Por otro 
lado, el FVS puede también autocali- 
brarse a nuevos datos si se incluyen atri- 
butos del arbolado, como índice de sitio, 
densidad y competencia. Las ecuaciones 
de crecimiento y las prescripciones de 
manejo también pueden ajustarse fácil- 
mente usando comandos específicos del 
FVS (Dixon, 2002). 

Para las simulaciones se usó un 
periodo de planeación de 100 años divi- 
dido en ciclos de 10 años. El periodo de 
100 años es el tiempo promedio que un 
árbol típico requiere para alcanzar la 
madurez en el área de estudio. La rege- 
neración de las especies dominantes 
Pinus arizonica y Quercos spp., se 
calculó en 840 y 346 árboles por hectárea 
a la edad de cinco años, respectiva- 
mente. El ciclo de corta para rodales inco- 
etáneos se determinó en 13 años, con 
una pendiente de la distribución diamé- 
trica igual a -1,4 (UAF Santiago Papas- 
quiaro, 2005). 

Usando el FVS se simularon nume- 
rosas condiciones y se examinaron los 
resultados al término del periodo de 
planeación (100 años). La solución 
óptima se determinó usando la ecuación 
4, substituyendo los valores resultantes 
de las simulaciones de cada una de las 
tres variables de decisión. Si cierta simu- 
lación daba valores consistentes para 
p= 2 (la distancia más corta entre dos 
puntos), los mismos valores de las varia- 
bles de decisión se calculaban al azar 
para p=1 y p ► <» para evaluar su 
consistencia. El criterio para determinar la 
mejor solución consistió en encontrar el 
mínimo valor de d p (la menor distancia al 
punto ideal) considerando la dinámica de 
crecimiento de los rodales en cada 
conglomerado. Para evaluar la variabi- 
lidad temporal de los conglomerados se 
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usaron mapas de distribución de estados 
de vegetación basados en el sistema de 
clasificación de dinámica de rodales de 
Oliver and Larson (1996). Las cuatro 
fases de vegetación, definidas parcial- 
mente por la magnitud de las variables de 
decisión, son etapa de iniciación del 
dosel, exclusión o eliminación por compe- 
tencia, reiniciación del dosel inferior y 
madurez. La etapa de madurez es la 
culminación del estado sucesional de la 
vegetación y comúnmente se le conoce 
como bosque clímax (ver detalles de esta 
clasificación en Oliver and Larson, 1996 y 
Mendoza et al., 2005). 

RESULTADOS Y DISCUSIÓN 

Los valores de la opinión de los 
expertos acerca de las condiciones 
iniciales del bosque con respecto a las 
tres variables de decisión ( ba , nth, y qmd ) 
se usaron para construir las funciones de 
respuesta para cada objetivo de manejo. 
Se usó una prueba de normalidad de los 
residuales para determinar la bondad de 
ajuste de los modelos, calculando el pará- 
metro Shapiro-Wilk y otros parámetros 
estadísticos como R 2 y nivel de signifi- 
cancia. Se calcularon las derivadas 
parciales de las diversas funciones obje- 
tivo para encontrar los valores máximos y 
mínimos. Los resultados del ajuste de 
modelos y valores máximos y mínimos se 
muestran en la tabla 2. En esta tabla se 
observa que los objetivos recreación, 
belleza escénica y productos no madera- 
bles tienen coeficientes negativos en 
cuanto a número de árboles por hectárea, 
lo cual sugiere que, en la opinión de los 
especialistas forestales consultados esta 
variable se ve favorecida en bosques 
abiertos. Esto coincide con los resultados 
de trabajos similares donde se ha repor- 
tado que, por lo general, la gente prefiere 
doseles abiertos y árboles grandes al 
momento de seleccionar sitios para la 
recreación (Tecle et al., 1998; Perez- 



Tabla 2. Ecuaciones de los objetivos, ajuste estadístico de modelos y valores locales mínimos y máximos^. 
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Verdín et al., 2008). Una situación similar 
fue observada para el objetivo venta de 
madera, donde bosques cerrados son 
generalmente percibidos como masas 
con árboles de diámetros menores, 
menos volumen y eventualmente menos 
ingresos por venta de productos madera- 
bles. 

Una vez construidas las funciones 
de respuesta y determinados los valores 
locales máximos y mínimos para todos 
los objetivos, se usó la ecuación 2 para 
determinar la solución compromiso para 
cada conglomerado en cada simulación. 
La tarea consistió en encontrar la menor 
distancia entre el juego de soluciones no 
dominadas y el punto ideal en cada simu- 
lación. Se usó un proceso iterativo para 
cada periodo de 10 años y para todos los 
conglomerados, simulando diferentes 
escenarios de manejo forestal hasta 
determinar la mejor solución. Además, se 
buscó obtener una mayor diversidad de 
estados de la vegetación al final del 
periodo de planeación. 

Este procedimiento no fue fácil. 
Algunos conglomerados se encontraban 
por abajo en términos de área basal, pero 
por arriba con respecto al número de 
árboles por hectárea. Además, el FVS no 
soporta más de una variable a la vez 
durante la simulación, así que en algunos 
conglomerados el proceso se repitió para 
las otras variables. Las simulaciones con 
el FVS incluyeron varios tratamientos 
silvícolas, como aclareos por lo bajo y por 
lo alto, selección individual, árboles 
padres, así como diferentes niveles de 
densidad residual. Tratamientos de no 
corta fueron desechados debido a la 
necesidad de seguir proveyendo un flujo 
constante de recursos económicos a los 
miembros del ejido. Esta diversidad de 
escenarios de manejo generó un esfuerzo 
considerable para estimar y comparar 
cada una de las d p . Los resultados de las 
simulaciones se analizaron mediante los 
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reportes estadísticos que incluyen, entre 
otros, los promedios de las variables de 
decisión y el sistema de visualización de 
rodales que genera esquemas tridimen- 
sionales de los resultados (Crookston y 
Dixon, 2005). 

La tabla 3 muestra las estructuras 
actuales y deseables para las tres varia- 
bles de decisión. Así por ejemplo, las 
estructuras iniciales del rodal tipo 1309 
del conglomerado uno son: 28,2 cm de 
diámetro cuadrático medio, 10,5 m 2 /ha 
de área basal, y 171 árboles por 
hectárea. Las estructuras deseables, de 
acuerdo a la información del inventario y 
la técnica de programación por compro- 
misos, para este mismo conglomerado 
son: 38,1 de diámetro cuadrático medio, 
17,2 m 2 /ha de área basal y 151 árboles 
por hectárea. En general, al aplicar 
promedios ponderados a todos los 
conglomerados, se observa que el área 
basal deseable es mayor que la actual 
en un 31%, el diámetro cuadrático medio 
es mayor en un 24% mientras que el 
número de árboles por hectárea es 
menor en un 5%. La tabla 3 muestra 
también el tratamiento inicial aplicado y 
los valores de la técnica de programa- 
ción por compromisos ajustados para 
p =2 (Ecuación 4). Estos valores corres- 
ponden a la distancia mínima entre la 
estructura actual y la deseable (punto 
ideal) para todos los objetivos mane- 
jados simultáneamente. La estructura 
deseable consideró no sólo una amplia 
variación de prescripciones silvícolas 
sino también la dinámica de crecimiento 
de los rodales. 

La figura 5 muestra la condición 
actual, a la mitad y al final del periodo de 
planeación (100 años) de las fases de 
vegetación. En la etapa inicial es evidente 
la dominancia de la fase de exclusión 
(88% del área forestal total) y en menor 
escala la fase de iniciación (12% del área 
forestal total). La fase de exclusión se 



Tabla 3. Estructuras deseadas y valores de la programación por compromisos (p= 2) en el Ejido Los Altares 
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Figura 5. Fases de la vegetación previstas a lo largo del periodo de planeación 
(mostradas para la condición actual, a la mitad y al final del periodo) 


caracteriza por una relativa dificultad para 
el establecimiento y desarrollo de nuevos 
individuos (Mendoza et al., 2005). A 
través de aclareos selectivos, el ejido y 
los servicios técnicos han mantenido esta 
fase a lo largo de los años para hacer 
más eficiente la producción maderable. A 
la mitad, y sobre todo al final del periodo 
de planeación, se pueden observar otras 
fases de la vegetación como reiniciación 
del dosel inferior y madurez. La presencia 
de otras fases vegetativas garantiza una 
mayor diversidad vegetal y animal, y 
como señalan Oliver y Larson (1996), es 
muy importante también para proveer 
valores recreacionales. Al final del 
periodo de planeación, los porcentajes de 
distribución de las fases vegetativas son: 
exclusión 51%, reiniciación del dosel infe- 
rior 24%, madurez con un estrato 11% y 
madurez con varios estratos 14%. 
Valores muy similares fueron encontrados 
por Mendoza et al. (2005) en una región 
del estado de Jalisco, donde la aplicación 


de varios tratamientos silvícolas ha gene- 
rado una mayor diversificación de fases 
vegetativas sin reducir el volumen de 
producción maderable. 

Análisis de sensibilidad 

Se realizó un análisis de sensibilidad 
para probar la fortaleza de los algoritmos 
con respecto a cambios en p y en el peso 
relativo de los objetivos. Dado que el 
proceso incluyó el examen simultáneo de 
tres variables de decisión continuas, se 
usó nuevamente la herramienta Solver 
del paquete Excel®. El primer paso de 
sensibilidad consistió en evaluar el algo- 
ritmo de la técnica de multicriterio con 
respecto a los valores de p= 1 , 2 e oc. Los 
resultados del análisis en este caso no 
mostraron diferencias en las soluciones 
excepto en los valores de nth cuando 
p = 1. Esto muestra que el algoritmo de 
multicriterio no es altamente sensible a 
cambios en los valores de p. 
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El segundo análisis se realizó para 
determinar la sensibilidad del algoritmo a 
cambios en los valores del nivel de prio- 
ridad (expresados en porcentaje) de los 
objetivos (Xt). En este caso, se separó la 
prioridad asignada por cada uno de los 
grupos participantes, es decir, los repre- 
sentantes del ejido Los Altares, los espe- 
cialistas forestales, funcionarios del 
gobierno y los integrantes de organiza- 
ciones no gubernamentales. De esta 
manera, el algoritmo fue evaluado con 
relación a la prioridad de los objetivos de 
manejo otorgada por cada grupo. Sin 
embargo, no hubo diferencias significa- 
tivas en las soluciones preferidas al modi- 
ficarse dichas prioridades. Los cuatro 
grupos participantes tienen preferencias 
similares en las ponderaciones asociadas 
con respecto a los objetivos del manejo, 
de tal forma que ninguno de ellos, indivi- 
dualmente o combinados, afectan los 
resultados. Para verificar que la solución 
no fue sensible a los cambios en las 
ponderaciones, se realizaron dos pruebas 
adicionales: (1) se incluyó un conjunto de 
prioridades asignadas de manera arbi- 
traria y (2) se aplicaron prioridades 
iguales para cada objetivo y luego se esti- 
maron las soluciones de compromiso bajo 
estos dos escenarios. En ambos casos, 
tampoco hubo cambios significativos en 
las soluciones compromiso. Por consi- 
guiente, los resultados obtenidos mues- 
tran que el algoritmo de la técnica de 
programación por compromisos no es 
sensible a cambios en los valores de p ni 
de Á,¡. Una de las razones por las que no 
se detectaron cambios sustanciales en el 
análisis de sensitividad puede deberse a 
las grandes diferencias entre los valores 
actuales y deseados de las estructuras 
forestales. Es decir, la distancia entre la 
frontera eficiente (dada por las estruc- 
turas actuales) y el punto ideal es tan 
amplia que cambios en p o X¡ no afectan 
significativamente las soluciones compro- 
miso. Si las estructuras forestales 
actuales, particularmente el área basal y 


el diámetro cuadrático medio fueran 
mayores, entonces las soluciones 
compromiso podrían verse afectadas por 
cambios en los valores de p y/o en la 
ponderación otorgada a los objetivos de 
manejo por los grupos participantes. 

RESUMEN Y CONCLUSIONES 

En este estudio se formuló y solu- 
cionó un problema de toma de decisiones 
aplicando técnicas multicriterio en el 
manejo integral de recursos naturales, 
usando como caso de estudio el ejido 
forestal Los Altares del estado de 
Durango. Con base en la opinión de 
especialistas forestales se evaluaron diez 
objetivos forestales y tres atributos daso- 
métricos. Para determinar el nivel de 
importancia (prioridades) de los objetivos 
de manejo, se consultaron a 144 
personas agrupadas en cuatro categorías 
que fueron: dueños y/o propietarios, 
especialistas forestales, representantes 
de gobierno y organizaciones no guber- 
namentales. Se realizaron numerosas 
simulaciones de manejo de los recursos 
forestales, evaluando el comportamiento 
dinámico de los tres atributos dasomé- 
tricos, con el fin de encontrar las mejores 
estructuras forestales que satisficieran los 
objetivos de manejo. Para determinar el 
juego de soluciones no dominadas, se 
usó la técnica de programación por 
compromisos que se ha utilizado con 
éxito para resolver problemas similares 
en el manejo de cuencas hidrográficas y 
recursos forestales. 

Los resultados de este trabajo 
indican que las estructuras forestales 
deseables, particularmente área basal y 
el diámetro cuadrático mínimo, general- 
mente estuvieron por encima de las 
condiciones actuales observadas en los 
rodales tipo. Es decir, los bosques 
actuales del ejido son en su mayoría 
rodales jóvenes con una alta densidad de 
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árboles por hectárea, de diámetros 
pequeños y área basal baja. Dichos 
bosques se encuentran mayoritariamente 
en una sola etapa vegetacional (exclu- 
sión) y carecen de rodales maduros y 
sobremaduros que puedan balancear la 
presencia de hábitat para la flora y fauna 
silvestre. Este tipo de estructuras difícil- 
mente puede facilitar el manejo de varios 
objetivos de manejo y proteger la inte- 
gridad del ecosistema. 

Muchos factores pueden contribuir a 
esas diferencias; por ejemplo, el manejo 
individual de un solo recurso (madera) 
puede reducir el nivel de existencias y la 
calidad de los árboles residuales gene- 
rando distribuciones diamétricas desequi- 
libradas (Gingrich, 1993). Eventualmente, 
este tipo de manejo forestal afecta la 
dinámica de crecimiento de las masas 
forestales al excluir la presencia de otras 
fases vegetativas. Excluyendo los incen- 
dios forestales que dentro del ejido se 
limitan a áreas pequeñas y de baja inten- 
sidad, otros factores que afectan la diná- 
mica de los bosques son la actividad 
humana, al deforestar áreas para la agri- 
cultura, caminos y pastoreo. La presencia 
de ganado, la cual es evidente en toda el 
área de estudio, frecuentemente afecta la 
regeneración. La interacción de todos 
estos factores genera riesgos para la inte- 
gridad y sustentabilidad de los ecosis- 
temas forestales. 

Los resultados sugieren que el 
manejo forestal enfocado a lograr obje- 
tivos múltiples es diferente al punto de 
vista tradicional, en el que la madera ha 
sido el producto más importante para 
maximizar ganancias. Pero lo más impor- 
tante es que el manejo forestal tradi- 
cional ha ignorado muchos valores no 
económicos, culturales y ambientales, 
como la biodiversidad, recreación y el 
mantenimiento del ciclo hidrológico. Se 
espera que la adopción del manejo de 
multirecursos y la dinámica de los 
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bosques genere una mayor colaboración 
entre los administradores del recurso 
forestal, los grupos ecologistas y otros 
interesados en el manejo de los ecosis- 
temas forestales para la producción de 
bienes y servicios ambientales. Se 
espera también que el ejemplo presen- 
tado aquí propicie un mejor entendi- 
miento, así como la aceptación y el uso 
de técnicas multicriterio en el manejo de 
los ecosistemas forestales. Las 
opiniones de expertos o consultores, si 
bien han probado su utilidad en muchos 
estudios, siguen siendo opiniones perso- 
nales y pueden reflejar aspiraciones, no 
necesariamente situaciones reales. Para 
mejorar la confiabilidad de este enfoque 
se requiere más investigación, especial- 
mente en la evaluación de los bienes y 
servicios ambientales que dependen de 
ciertos atributos del rodal. Es necesario 
el establecimiento de parcelas experi- 
mentales para generar más conocimiento 
en las variables para evaluar los bienes y 
servicios que ofrece el recurso forestal. 
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Variación en la cobertura de suelo 
en un ensayo de procedencias de 
Pinus greggii Engelm. 
en el cerro El Potosí, Galeana, Nuevo León 

Variations on soil protection in a test of provenances of 
Pinus greggii Engelm. in El Potosí Hill, Galeana, Nuevo León 

Rodrigo Rodríguez-Laguna 1 , Joel Meza-Rangel 1 , 

Jesús Vargas-Hernández 2 y Javier Jiménez-Pérez 3 

RESUMEN 


Con la finalidad de conocer la protección al suelo que brindan las copas de los árboles en un 
ensayo de procedencias de Pinus greggii Engelm. a 4.5 años de plantados en el Cerro El Potosí, 
Galeana, Nuevo León, se midió supervivencia, altura y características de la copa; también se calcu- 
laron área de proyección de copa, área de intercepción lumínica y porcentaje de copa. No se encon- 
traron diferencias significativas entre procedencias en la supervivencia, pero se tuvo en general 93% 
de supervivencia; al contrarío sí hubo diferencias significativas (P<0,05) en altura, altura a la primera 
rama, altura al punto amplio de la copa, área de proyección de copa, área de intercepción lumínica y 
porcentaje de copa. Los árboles de las procedencias con mayor área de proyección de copa tuvieron 
0,85 m 2 y en el área de Intercepción lumínica varió en más de 1 ,1 m 2 entre las procedencias extremas. 
Para el porcentaje de copa el mayor valor obtenido fue de 89,1% de tallos con copa. Los componentes 
de varianza mostraron que las procedencias aportaron 7,5% de la variación total en el área de proyec- 
ción de copa y en el diámetro de copa. En cambio, este valor disminuyó para la altura a la primera 
rama y la altura al punto amplio de la copa con 5,2% y 2,9%, respectivamente. A 4.5 años de estable- 
cido el ensayo la copa de los árboles ha cubierto el 14,8% de la superficie total plantada. El creci- 
miento de los árboles se correlacionó positivamente con el área de intercepción lumínica (r = 0,94), 
área de proyección de copa (r = 0,86) y porcentaje de copa (r = 0,48). El estudio muestra evidencias 
para seleccionar la procedencia de Los Lirios, Coahuila y Agua Fría, Nuevo León, para proteger el 
suelo del sitio de plantación. 

PALABRAS CLAVE: 

Cobertura de copas, Pinus greggii Engelm., supervivencia. 


ABSTRACT 

With the purpose of knowing the soil protection that tree crowns provide in a test of provenances 
of Pinus greggii Engelm. planted 4.5 years prior to the test in the Cerro El Potosí, Galeana, Nuevo 
León, survival, heíght and crown characteristics were measured; crown projection area, luminance 
interception area and percentage of crown were also calculated. No significant difference was found 
between provenances in the survival, but in general it was found a survival rate of 93%; however, there 
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were significant differences (P<0,05) in 
height, height at the first branch, height 
to the wide point of the crown, crown 
projection area, luminance ¡nterception 
area and percentage of crown. The trees 
of the provenances with greater crown 
projection area had 0,85 m 2 and in the 
luminance ¡nterception area varied in 
more of 1,1 m 2 between the prove- 
nances extremes. For the crown percen- 
tage the higher valué obtained was of 
89,1% of stems with crown. The compo- 
nents of variance showed that prove- 
nances provided 7.5% of the total varia- 
tion in crown projection area and in the 
crown diameter. On the other hand, this 
valué diminished for height at the first 
branch and height to the wide point of 
the crown with 5,2% and 2,9%, respecti- 
vely. After 4,5 years of establishing the 
test, the crowns of trees have covered 
14,8% of the total surface planted. The 
tree growth was positively correlated 
with the luminance ¡nterception area (r = 
0,94), crown projection area (r = 0,86) 
and percentage of crown (r = 0,48). The 
study shows evidences to select the 
provenance of Los Lirios, Coahuila and 
Agua Fría, Nuevo León as the best 
choices to protect the soil at the planting 
site. 

KEY WORDS: 

Crowns coverage, Pinus greggii Engelm., 
survival. 


INTRODUCCIÓN 

En México existe la tendencia a 
transformar las mejores áreas de voca- 
ción forestal hacia áreas de uso agrícola 
y pecuario. Lamentablemente dichas 
áreas en pocos años pierden su fertilidad 
y alcanzan rendimientos cada vez 
menores. La degradación de áreas en los 
diversos ecosistemas forestales trae 
como consecuencia la reducción de la 
biodiversidad, lo cual se refleja en la 
disminución de bienes y servicios que 
brindan a la sociedad (Hudson, 1982). 
Los ecosistemas forestales son afectados 
por diversos factores que inciden en su 
degradación y deforestación, los princi- 
pales son los desmontes con 94%, incen- 
dios forestales con 2.1%, el cambio de 


uso del suelo con autorización con 1 .3%, 
las talas clandestinas con 1.0%, plagas y 
enfermedades forestales con 0.5% y 
otros con 1.0% (SEMARNAT, 2000). A lo 
anterior hay que añadir la acción del 
hombre en los riesgos de degradación, 
debido a los efectos de las perturba- 
ciones asociadas con el manejo intensivo 
que se hace en los suelos sin un conoci- 
miento adecuado de su dinámica, 
llevando consigo un deterioro acelerado 
del ecosistema (González et al., 1990). 

Los ensayos de procedencias por 
especie es el procedimiento experimental 
más empleado para comparar el compor- 
tamiento y la productividad de semilla que 
proviene de diferentes poblaciones. Estos 
ensayos se pueden establecer en 
ambientes extremos en los cuales se 
pretende establecer las plantaciones, con 
el argumento de encontrar las mejores 
fuentes de semilla, en términos de adap- 
tación al sitio de plantación, conformación 
de las copas y en la productividad de los 
árboles (Zobel y Talbert, 1992). 

Pinus greggii Engelm. es una 
especie forestal nativa de México, que se 
distribuye en pequeños rodales a lo largo 
de la Sierra Madre Oriental (Martínez, 
1992). Esta especie se considera rústica 
y con buenas características para 
desarrollarse sobre terrenos pobres 
(González, 1978), posee la habilidad de 
crecer rápidamente y ha sido reportado 
como tolerante tanto a la sequía (Vargas 
y Muñoz, 1988) como al ataque de 
algunas plagas fuera de su área natural 
(Dvorak et al., 2000). De igual forma 
estas características le han hecho una 
especie muy utilizada con fines de protec- 
ción o recuperación de áreas degradadas 
en las últimas décadas; algunos estudios 
demuestran que puede establecerse 
fuera de su área de distribución natural 
con resultados favorables (Dvorak et al., 
1996; Velasco, 2001; Valencia et al., 
2006). 
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La capacidad de protección o recu- 
peración de un suelo depende de la velo- 
cidad con que se logre la cobertura de 
éste por los árboles; ya que al aumentar 
la cobertura de copa se reduce el impacto 
de la lluvia y de otros agentes erosivos 
sobre el suelo, se aumenta captación de 
energía y la productividad primaria, así 
como la aportación de materia orgánica al 
suelo. Estos procesos ayudan a formar 
un microclima propicio para los agentes 
bióticos que incorporan la materia orgá- 
nica al suelo, dejándola en condiciones 
para ser utilizada nuevamente por las 
plantas (Tang et al., 1999). 

La copa del árbol es el componente 
de la producción primaria y sus dimen- 
siones reflejan la salud del individuo; por 
ejemplo, las copas densas y largas están 
asociadas con un crecimiento vigoroso; 
las copas de bajo desarrollo y poco 
densas se presentan en sitios en condi- 
ciones desfavorables por la competencia, 
estrés por humedad o la influencia de la 
defoliación por insectos y enfermedades 
de las hojas (Schomaker et al., 1999). 
Algunos estudios muestran la importancia 
de conocer las características de la copa 
de los árboles y con ellas predecir la 
respuesta en crecimiento de los árboles 
(Doruska y Burkhart, 1994; Brunner, 

1998) , del mismo modo autores 
mencionan que la cantidad de luz que 
intercepta la copa determina en gran 
medida el crecimiento del árbol (Aiba y 
Kohyama, 1997; Sterck, 1999; Tang etal., 

1999) . 

Contrariamente a la importancia de 
proteger y evitar la erosión en el suelo, 
existe escasa información referente a la 
cobertura de suelo por las copas de los 
árboles en estados juveniles, en sitios 
con posibilidad de restauración. Si se 
conoce el tamaño y la forma de la copa de 
los árboles de las distintas procedencias 
es posible seleccionar árboles de esta 
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especie que cubran mayor superficie en 
menos tiempo. 

Tomando en cuenta los antece- 
dentes señalados se planteó como primer 
objetivo determinar la variación por proce- 
dencia en la supervivencia, cobertura de 
copa, superficie de intercepción lumínica 
y tamaño de copa, además de calcular la 
superficie de suelo cubierta por las copas 
de los árboles en un ensayo de nueve 
procedencias de Pinus greggii Engelm. a 
4.5 años de plantados en el cerro El 
Potosí, Galeana, Nuevo León. El segundo 
objetivo fue encontrar la relación entre 
cobertura de la copa, área de intercep- 
ción lumínica y tamaño de la copa con el 
crecimiento de los árboles. 


METODOLOGÍA 

Origen del germoplasma y ubicación 
del ensayo 

El estudio se realizó en un ensayo 
de P. greggii de 4.5 años de plantados en 
campo, en donde se evalúa germoplasma 
de 9 poblaciones naturales de P. greggii 
de los estados de Coahuila y Nuevo 
León. La ubicación geográfica y los datos 
más importantes de los sitios de colecta 
(procedencias) se presentan en la tabla 1 . 
La semilla se colectó en 1997 y se 
sembró en noviembre de 1998, en 
envases individuales de plástico negro de 
125 mi de capacidad dentro de un inver- 
nadero en el Colegio de Postgraduados. 
Se utilizó sustrato con 60% peat moss, 
30% vermiculita, 10% agrolita y 4 kilo- 
gramos de Osmocote 17-17-17 por metro 
cúbico de sustrato. Las plantas se mantu- 
vieron con riegos programados (3 riegos 
por semana); a los tres meses de edad se 
redujo a la mitad la densidad de plantas y 
se sacaron a platabandas con malla 
sombra. 
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Tabla 1. Información de las procedencias de la semilla de P. greggii utilizadas. 


Procedencia 

Coordenadas geográficas 
Lat. Norte Long. Oeste 

Altitud 

(msnm) 

Precipitación 

(mm) 

Puerto Ei Conejo, N. L. 

25° 28’ 

100° 35’ 

2520 

650 

Santa Anita, Coah. 

25° 27’ 

100° 34’ 

2560 

650 

Agua Fría, N. L. 

25° 26’ 

100° 30’ 

2400 

633 

Puerto San Juan, Coah. 

25° 25’ 

100° 33’ 

2613 

600 

Los Lirios, Coah. 

25° 23’ 

100° 31’ 

2420 

600 

El Penitente, Coah. 

25° 22’ 

100° 54’ 

2405 

500 

Jamé, Coah. 

25° 21’ 

100° 34’ 

2552 

600 

Las Placetas, N. L. 

24° 55’ 

100° ir 

2450 

750 

La Tapona, N. L. 

24° 43’ 

100° 06’ 

2130 

650 


El ensayo de procedencias se esta- 
bleció en septiembre de 1999 en una 
parcela del Ejido 18 de Marzo, Galeana, 
Nuevo León. En los 24° 53’ N y 100° 12’ 
O, a una altitud de 2209 msnm; la precipi- 
tación media anual es de 401.1 mm y la 
temperatura media anual es de 18.2 °C, 
se considera un clima seco semi-cálido, 
invierno fresco, muy extremoso, con 
lluvias de verano y precipitación invernal 
no muy significativa. Los meses más 
lluviosos en la región son agosto y 
septiembre (García, 1988). 

Condiciones del terreno en el ensayo 

El ensayo se estableció en un 
terreno agrícola abandonado, con una 
pendiente mayor de 19%, con exposición 
noreste. Debido al uso agrícola anterior, 
el suelo presentaba un nivel de erosión 
moderado. Para evitar una erosión mayor 
se hicieron cepas de 30x30x30 cm, esta- 
bleciendo la planta en el centro de la cepa 
y formando un cajete para captar agua de 


lluvia. La textura del suelo en la capa 
superficial es limosa-arcillosa, con estruc- 
tura terregosa moderadamente dura. 

Diseño experimental del ensayo 

Para controlar el componente 
ambiental, especialmente las condiciones 
de suelo y de pendiente, en el ensayo se 
utilizó un diseño de bloques completos al 
azar, con 30 bloques transversales a la 
pendiente. Las nueve procedencias se asig- 
naron al azar dentro de cada bloque, las 
parcelas quedaron en hilera de cuatro 
plantas a lo largo de la pendiente. Las 
plantas se establecieron a un espaciamiento 
de 2x2 metros para generar una densidad 
de 2500 plantas por hectárea, ocupando el 
ensayo una superficie total de 4320 m 2 . 

Variables medidas 

En los árboles después de 4.5 años 
de plantados en campo, se determinó el 
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porcentaje de supervivencia en cada 
parcela. En esa misma fecha se midió 
altura total, altura a la primera rama viva, 
altura al punto amplio de la copa y el 
diámetro de la copa en direcciones norte- 
sur y este-oeste. 

Con los datos se procedió a calcular 
otras características de la copa del árbol, 
como: porcentaje de copa (PC), área de 
proyección de copa en el suelo (APCS) y 
área de intercepción lumínica de la copa 
(AILC) (Jiménez et al., 2002). El PC se 
define como la relación entre la longitud 
de copa (LC) y la altura total del árbol (Ht) 
expresado en porcentaje o como una 
fracción (Figura 1). Para determinar el 
porcentaje de copa se utilizó la ecuación 
siguiente: 


r r 

PC = —x 100 
Ht 


El APCS se estimó como la proyec- 
ción horizontal de la copa del árbol en el 
suelo, y se utilizó como una medida de 
densidad de la plantación (Figura 1): 


APCS = —*DPC 2 

4 


Donde DPC es el diámetro promedio 
de la copa del árbol. 


El AILC se conoce como la región 
fotosintética más activa de la copa en 
donde se absorbe la radiación lumínica y 
comprende la superficie lateral de la copa 
eliminando la superficie de sombra de la 
copa (Jiménez et al., 2002), se expresa 
en m 2 (Figura 1): 


AILC = 


jt *r 
6LC 2 



í 



Donde LC = longitud de copa del 
árbol (m) y r = es el radio de la copa (m) 


Altura T otal 



Figura 1. Variables medidas en campo y variables calculadas por árbol. 
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Análisis estadístico 

Se sometieron a un análisis de 
varianza con los promedios por parcela 
las variables: supervivencia, altura a la 
primera rama, altura al punto amplio de la 
copa, así como las variables calculadas 
de la copa. En el caso de la supervi- 
vencia, previo al análisis de varianza los 
datos fueron transformados con la función 
arco-seno, pero los valores promedio de 
las procedencias fueron expresados en 
las unidades originales. A las variables 
que presentaron diferencias estadísticas 
con P < 0.05 en el análisis de varianza se 
hizo la prueba de comparación múltiple 
de medias de Tukey. Con los resultados 
del análisis de varianza se estimó la 
contribución de las procedencias a la 
varianza total observada en cada 
variable, a partir de la estructura de la 
esperanza de los cuadrados medios del 
modelo estadístico para cada caso. 
También con los valores promedio por 
procedencia se estimó la correlación 
simple entre las variables evaluadas de la 
copa con el crecimiento en altura de los 
árboles, para conocer las posibles rela- 
ciones en la plantación. 

RESULTADOS Y DISCUSIÓN 
Supervivencia 

El análisis de varianza mostró que 
no hubo diferencias estadísticas significa- 
tivas (P < 0.05) entre procedencias en la 
supervivencia de las plantas a 4.5 años 
de haber establecido el ensayo; en 
general el ensayo tuvo una supervivencia 
de 92.7%, con valores promedio 
extremos de 96.4% para Puerto San Juan 
y de 86% para Agua Fría. 

En otros ensayos de P. greggii esta- 
blecidos en regiones con condiciones 
ambientales muy diversas tampoco se 
han encontrado diferencias significativas 


entre procedencias en la supervivencia 
de las plantas (Dvorak et al., 1996; 
Órnelas et al., 2001; Velasco, 2001; 
López-Upton et al., 2004; Contreras, 
2005), lo que muestra la plasticidad de la 
especie en el centro y norte del país. 
Según Zobel y Talbert (1992) las variables 
de adaptación como supervivencia, 
tienen un fuerte efecto de procedencia u 
origen geográfico, pero esto no se reflejó 
en el presente estudio. Es posible que el 
grupo de procedencias de P. greggii en el 
norte del país haya evolucionado de 
manera semejante, al estar expuesto a 
condiciones ambientales similares. 

El ensayo tuvo un grupo de seis 
procedencias con promedios por arriba 
de 93% de supervivencia (Puerto San 
Juan, Las Placetas, Puerto El Conejo, 
Los Lirios, La Tapona y Jamé) como se 
muestra en la figura 2. El valor promedio 
en otros estudios con la misma especie y 
utilizando algunas de las mismas proce- 
dencias tuvieron resultados menos satis- 
factorios como el ensayo realizado en 
Patoltecoya, Puebla; la supervivencia 
promedio fue de 75% a los dos años siete 
meses (López-Ayala et al., 1999) y de 
60% en el mismo ensayo a los seis años 
(López-Upton et al., 2004); en otro 
ensayo de procedencias realizado en el 
CAESA, Arteaga, Coahuila, se encontró 
68.6% de supervivencia a los cinco años 
con siete meses (Órnelas et al., 2001). 

Altura a la primera rama y al punto 
amplio de la copa 

En características de la copa como 
altura a la primera rama y altura al punto 
amplio de la copa se encontraron diferen- 
cias significativas (P < 0.05) entre las proce- 
dencias (Tabla 2). La característica que 
define la longitud de la copa de un árbol, es 
la altura de la primera rama de la copa, en 
este caso Puerto San Juan tiene la primer 
rama de la copa a ras de suelo con un valor 
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Figura 2. Supervivencia promedio de nueve procedencias de P greggii a 4.5 años de 
plantadas. Las barras con líneas diagonales muestran las procedencias de Coahuila, 

las barras con puntos son de Nuevo León. 


Tabla 2. Resultados del análisis de varianza y significancia en variables de la copa del 

árbol en nueve procedencias de P. greggii. 




Cuadrados Medios 


Variables 

Bbques 

Procedencia 

Bloq*Proc 

Error 


(29) a 

(8) 

(231) 

(651) 

Altura a la primera rama 

177.3 

316.8** 

112.6 

95.7 

Altura al punto amplio de la copa 

3015.1 

922.5* 

495.5 

240.3 


** 


En paréntesis se presentan los grados de libertad correspondientes a cada fuente de variación. 
Significativo al P<0.05 
Significativo al P<0.01. 
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de 1 1 .4 cm y la menor altura al punto amplio 
de la copa (54.2 cm), disminuyendo el 
impacto de la gotas de lluvia en el suelo y 
formando un microclima que le ayuda a 
mantener la humedad por mayor tiempo, en 
contraste Los Lirios tuvo la primera rama a 
17.5 cm del suelo y Agua Fría tuvo 10.1 cm 
más arriba el punto amplio de la copa (64.3 
cm) (Figura 3). Si bien el sitio de plantación 
se encuentra dentro del área de distribución 
natural de la especie, se han registrado dife- 
rencias entre procedencias en caracterís- 
ticas de la copa, en respuesta a las condi- 
ciones de humedad, temperatura y suelo. 

Los componentes de varianza en la 
altura a la primera rama muestran que el 
5.2% de la variación total se debe a las 
procedencias; de la misma forma, la 
variación total en la altura al punto amplio 
de la copa es de 2.9% que se debe a las 
procedencias. El mayor porcentaje de 
variación está explicado por el error 
estándar, que se atribuye a factores que 
no pudieron ser controlados en el ensayo. 

Diámetro y porcentaje de copa 

Las procedencias de P. greggii 
difieren en forma significativa (P < 0.01) 


en el diámetro de copa y porcentaje de 
copa de los árboles (Tabla 3). La proce- 
dencia de Agua Fría tiene el mayor 
diámetro promedio de copa, con 99.2 cm, 
mientras que la procedencia de Puerto 
San Juan se ubica en el extremo inferior 
con 79.9 cm. Para el porcentaje de copa, 
los árboles de Agua Fría tuvieron 89.1% 
de tallos con copa, en tanto que Santa 
Anita tuvo el valor extremo inferior con 
82.4% (Figura 4). Puede ser que las 
procedencias que en este momento 
presentaron mayor porcentaje de copa, 
después tengan una reducción significa- 
tiva debido a la autopoda que presenta la 
especie en condiciones de competencia 
por luz y espacio. Los componentes de 
varianza muestran que 7.5 % de la varia- 
ción total en el diámetro de copa del árbol 
se debe a las procedencias; en cambio, 
para el porcentaje de copa del árbol este 
valor se redujo a 4.4%. 

Valencia etal. (2006), encontraron que 
el diámetro de la copa varió de 46.9 a 74.9 
cm en árboles de Pinus greggii de 2.5 años 
de edad en Tlacotepec Plumas, y de 43 a 
57.5 cm en Magdalena Zahuatlán, en la 
Mixteca Alta de Oaxaca. Hay escasa infor- 
mación sobre las características de la copa 



Figura 3. Valores promedio en (A) altura a la primera rama; (B) altura al punto amplio de 
la copa en nueve procedencias de P. greggii a 4.5 años de plantadas. Barras con la 

misma letra son estadísticamente iguales. 
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Tabla 3. Resultados del análisis de varianza para diámetro de copa y porcentaje de 

copa en nueve procedencias de P. greggii 




Cuadrados Medios 


Variables 





Bloques 

Procedencias 

Bloq*Proc 

Error 


(29) a 

(8) 

(231) 

(651) 

Diámetro de copa 

6425.2 

3845.1** 

1127.1 

548.7 

Porcentaje de copa 

338.1 

350.4* 

154.5 

131.1 


En paréntesis se presentan los grados de libertad correspondientes a cada fuente de variación. 
Significativo al P<0.05 
Significativo al P<0.01. 



Figura 4. Valores promedio en (A) diámetro de copa; (B) pocentaje, de la copa en nueve 
procedencias de P. greggii a 4.5 años de plantadas. Barras con la misma letra 

son estadísticamente iguales. 


en edades tempranas de esta especie. En 
este ensayo del cerro El Potosí, el diámetro 
de copa varió de 79.9 a 99.2 cm. 

Área de proyección de copa y área de 
intercepción lumínica 

El análisis de varianza mostró dife- 
rencias significativas (P < 0.01) entre las 
procedencias en el área de proyección de 
copa y en el área de intercepción lumínica 
(P < 0.10) (Tabla 4). En promedio, cada 


árbol de Los Lirios cubrió 0.85 m 2 de 
terreno con su copa, mientras que los de 
Puerto San Juan sólo cubrieron 0.55 m 2 
(Tabla 5), lo que representa 35.3% menos 
de cobertura por árbol. Esta diferencia en 
los primeros años de vida en referencia a 
la cobertura de suelo implica tiempo para 
recuperar la cobertura total de la vegeta- 
ción. Las procedencias con mayor cober- 
tura permitirían disminuir más pronto el 
problema de la erosión de suelo, además 
de aportar más residuos vegetales y favo- 
recer el mejoramiento del suelo. 
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Tabla 4. Análisis de varianza para el área de proyección de copa y el área de 
intercepción lumínica en nueve procedencias de P. greggii. 


\/ar¡ahlp“; 


Cuadrados Medios 



Bloques 

(29) a 

Procedencias 

(8) 

Bloq*Proc 

(231) 

Error 

(651) 

Área de proyección de copa 
sobre el suelo 

1,47 

0,98** 

0,27 

0,13 

Área de intercepción 
lumínica de la copa 

25,4 

9,5* 

5,1 

1,93 


En paréntesis se presentan los grados de libertad correspondientes a cada fuente de variación. 
Significativo al P<0.10 
Significativo al P<0.01. 


Los componentes de varianza 
mostraron que 7.5% de la variación total 
en la cobertura de la copa se debe a las 
procedencias, mientras que 92.5% de la 
variación se encuentra dentro de las 
parcelas chicas (error de muestreo). Lo 
anterior indica que existe una amplia varia- 
ción en la cobertura de la copa sobre el 
suelo a nivel de procedencias, lo cual es 
de gran interés en la selección de material 
en un programa de mejoramiento genético 
con la finalidad de protección al suelo. 

Los árboles de Agua Fría presen- 
taron la mayor área de intercepción lumí- 
nica promedio (3.0 m 2 ) que superó a los 
de Puerto San Juan (1.9 m 2 ). Los árboles 
con mayor área de intercepción lumínica 
podrían tener mayor potencial de acti- 
vidad fotosintética y, por tanto, mayor 
potencial de crecimiento. Schomaker et 
al. (1999), indicaron que la copa es uno 
de los principales componentes de la 
producción primaria en los árboles y sus 
dimensiones reflejan el vigor de los indivi- 
duos. Cuando la finalidad de la plantación 
es dar protección al suelo, es de utilidad 
identificar las procedencias que desarro- 
llan copas más amplias con mayor 
rapidez. Los componentes de varianza 
mostraron que sólo 2.7% de la variación 
total del área de intercepción lumínica se 
debe a las procedencias. 


Después de 4.5 años de plantados 
los árboles en conjunto de las nueve 
procedencias tuvieron una cobertura de 
suelo de 638.1 m 2 , lo que representa 
14.8% del total de la superficie. La proce- 
dencia de Los Lirios aporta 13.9% de la 
cobertura total del suelo en el sitio, en 
tanto que Puerto San Juan sólo cubre 
9.3% (Tabla 5). Lo anterior indica que 
existe un potencial para seleccionar 
procedencias que protejan mayor super- 
ficie de suelo desde los primeros años. 

Análisis de correlación 

El crecimiento de los árboles se corre- 
lacionó positivamente (P < 0.01) con el 
área de intercepción lumínica de la copa 
(r = 0.94**), con el área de proyección de 
copa (r = 0.86**) y con el porcentaje de 
copa del árbol (r = 0.48**). Lo anterior signi- 
fica que las procedencias con mayor cober- 
tura del suelo, tuvieron árboles con mayor 
altura y mayor cantidad de área de inter- 
cepción lumínica; en consecuencia, los 
árboles presentaron un mayor porcentaje 
de copa. Lo anterior apoya lo que 
mencionan González (1978), Eguiluz 
(1978) y Perry (1991) al establecer que P. 
greggii es una especie muy rústica y que 
posee buenas características para desarro- 
llarse con cierta facilidad sobre terrenos 
pobres y presentar un rápido crecimiento. 
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Tabla 5. Valores promedio por procedencia (m 2 , %) de cobertura del suelo en nueve 
procedencias de P greggii a 4.5 años de plantadas. 


Procedencia 

Área de proyección de copa 
sobre el suelo por árbol (m 2 ) 

Área de proyección de 
copa sobre el suelo por 
procedencia (m 2 ) 

Área de proyección de 
copa sobre el suelo 

(%) 

Puerto El Conejo, N. L. 

0,59 

61,4 

9,6 

Santa Anita, Coah. 

0,69 

67,6 

10,6 

Agua Fría, N. L. 

0,81 

77,0 

12,1 

Puerto San Juan, Coah. 

0,55 

59,4 

9,3 

Los Lirios, Coah. 

0,85 

88,4 

13,9 

El Penitente, Coah. 

0,74 

71,0 

11,1 

Jamé, Coah. 

0,66 

68,0 

10,7 

Las Placetas, N. L. 

0,68 

72,8 

11,4 

La Tapona, N. L. 

0,69 

72,5 

11,4 

Total 


638,1 m z 

100 % 


CONCLUSIONES 

Existe variación significativa entre 
procedencias de Pinus greggii Engelm. 
en el área de proyección de copa, en el 
área de intercepción lumínica y en el 
porcentaje de copa de los árboles. Pero 
no difirieron en la supervivencia y obtu- 
vieron en general 93% de supervivencia 
en las procedencias. 

Las nueve procedencias de Pinus 
greggii Engelm. del norte del país 
cubrieron 14.8% de la superficie del suelo 
plantada. No obstante, las procedencias 
que presentaron mayores ventajas para 
usarse en plantaciones con protección al 
suelo en la región sur de Nuevo León, 
fueron Los Lirios, Coahuila y Agua Fría, 
Nuevo León, porque sus árboles obtu- 
vieron mayores valores en características 
de la copa a los 4.5 años de plantados. 

Las procedencias que tienen árboles 
con mayor área de proyección de copa, 
área de intercepción lumínica y porcen- 


taje de la copa tuvieron correlación posi- 
tiva con el crecimiento de los árboles, lo 
que significa que las procedencias con 
mayor cobertura del suelo, tuvieron 
árboles con mayor altura y también mayor 
cantidad de área de intercepción lumínica 
y con copas más grandes. 

AGRADECIMIENTOS 

Al Fondo Mexicano para la Conser- 
vación de la Naturaleza, por aportar los 
recursos financieros para la colecta de 
semilla y el establecimiento de la planta- 
ción dentro del proyecto de investigación 
“Diversidad genética, demografía y 
conservación de poblaciones naturales 
de Pinus greggii Engelm.”, dirigido por el 
Dr. Jesús Vargas Hernández. También se 
agradece ampliamente a los revisores de 
la revista Madera y Bosques que permi- 
tieron mejorar el manuscrito. 


58 


Variación en la cobertura de suelo en un ensayo de procedencias de Pinus greggii 


REFERENCIAS 

Aiba, S. y T. Kohyama. 1997. Crown 
architecture and life-history traits at 
14 tree species in a warm tempera- 
ture rain forest significance of spatial 
heterogeneity. Journal of Ecology 85 
(5): 611-624. 

Brunner, A. 1998. A light model for 
spatially explicit forest stands 
models. Forest Ecology and Mana- 
gement 147 (1-3): 19-46. 

Contreras, M. R. 2005. Ensayo de tres 
procedencias de Pinus greggii 
Engelm. establecido en el C. A. E. S. 
A. Arteaga, Coahuila. Tesis profe- 
sional. UAAAN. Buenavista, Saltillo, 
Coahuila. 58 p. 

Doruska, D.F. y H.E. Burkhart. 1994. 
Modeling the diameter and locational 
distribution of branches within the 
crowns of loblolly pine trees in 
unthinned plantations. Canadian 
Journal or Forestry Research 24: 
2362-2376. 

Dvorak, W.S., J.E. Kietzka y J.K. 
Donahue. 1996. Three-year survival 
and growth of pro ve nances of Pinus 
greggii in the tropics and subtropics. 
Forest Ecology and Management 83: 
123-131. 

Dvorak, W.S., J.E. Kietzka, J.K. Donahue, 
G.R. Hodge y T.K. Stanger. 2000. 
Pinus greggii, In: Conservaron & 
Testing of Tropical & Subtropical 
Forest Tree Species by the 
CAMCORE Cooperative. College of 
Natural Resources, NCSU. Raleigh. 
NC. USA. p. 52-73. 

Eguiluz, R T. 1978. Ensayo de integración 
de los conocimientos del género Pinus 
en México. Tesis de Licenciatura. 
Universidad Autónoma Chapingo. 
Chapingo, México, p. 438-446 


García, E. 1988. Modificaciones al 
sistema de clasificación climática de 
Kóppen. Cuarta edición. México, 
D.F. 219 p. 

González, V.,C.E. 1978. Breve análisis de 
la investigación sobre plantaciones 
forestales de la Dirección General 
de Investigación y Capacitación 
Forestal. INIFAP. Publicación espe- 
cial No. 13. p. 83-87. 

González, K.V., G. Hernández S. y D. 
Hernández S. 1990. Evaluación de 
la pérdida de suelo (Prmera aproxi- 
mación) en una región de la serranía 
de Actopan, Hidalgo, México. In: 
Memorias del Primer Simposium 
Nacional “Degradación del suelo”. 

Hudson, N. 1982. Conservación del 
suelo. Ed. Reverté. Barcelona, 
España. 335 p. 

Jiménez, J., H. Kramer y O. Aguirre. 
2002. Bestandesuntersuchungen in 
einem ungleichaltrigen Tannen-, 
Douglasien-, Kiefern- Naturbestand 
Nordostmexikos. Allgemeine Forst 
und Jagdzeitung 173: 47-55. 

López-Ayala, J.L., J.J. Vargas H., C. 
Ramírez H. y J. López U. 1999. 
Variación intraespecífica en el patrón 
de crecimiento del brote terminal de 
Pinus greggii Engelm. Revista 
Chapingo serie Ciencias Forestales 
y del Ambiente 5(2): 133-140. 

López-Upton, J., C. Ramírez H., O. Plas- 
cencia E. y J. Jasso M. 2004. Variación 
en crecimiento de diferentes pobla- 
ciones de las dos variedades de Pinus 
greggii. Agrociencia 38(4): 457-464. 

Martínez, M. 1992. Los pinos mexicanos. 
Tercera edición. Ed. Botas. 368 p. 

Órnelas, H.G., E. Aldrete M. y E.H. 
Cornejo O. 2001. Ensayo de tres 


Madera y Bosques 15(1), 2009:47-59 

procedencias de Pinus greggii 
Engelm. en el Campo Agrícola Expe- 
rimental Sierra de Arteaga. Foresta- 
AN. Nota Técnica No. 6. UAAAN. 
Saltillo, Coah. 12 p. 

Perry, J.P. Jr. 1991. The pines of México 
and Central America. Timber Press. 
Portland, Oregon. 231 p. 

Schomaker, M., S. Zanoch y K. Stolte. 

1999. Tree crown condition indicator. 
USDA. Forest Service. Forest Health 
Monitoring Fact Sheet. 

SEMARNAT (Secretaría del Medio 
Ambiente y Recursos Naturales). 

2000. El sector forestal de México: 
Situación actual y perspectivas 
(1998). Sistema Nacional de Infor- 
mación Forestal. 

Sterck, F.J. 1999. Crown development in 
tropical rain forest trees in gaps and 
understorey. Plant Ecology 143: 89- 
98. 

Tang, Z., J.L. Chambers, S. Duddanti, S. Yu 
y J.P. Barnett 1999. Seasonal shoot 
and needle growth of loblolly pine 
responds to thinning, fertilization and 
crown position. Forest Ecology and 
Management 120: 117-130. 


59 

Valencia, M.S., M.V. Velasco G., M. Gómez 
C., M. Ruiz M. y M.A. Capó A. 2006. 
Ensayo de procedencias de Pinus 
greggii Engelm. en dos localidades de 
la Mixteca Alta de Oaxaca, México. 
Rev. Fitotec. Mex. 29: 27-32. 

Vargas, H.J.J. y A. Muñoz O. 1988. Resis- 
tencia a sequía II. Crecimiento y 
supervivencia en plántulas de cuatro 
especies de Pinus. Agrociencia 72: 
197-208. 

Velasco, G.M.V. 2001. Ensayo de 13 
procedencias de Pinus greggii 
Engelm. en dos localidades de la 
Mixteca Alta, Oaxaca. Tesis profe- 
sional. UAAAN. Buenavista, Saltillo, 
Coahuila, México. 75 p. 

Zobel, B. y J. Talbert. 1992. Técnicas de 
mejoramiento de árboles forestales. 
Versión española por Manuel 
Guzmán O. Primera reimpresión. 
Editorial LIMUSA. México. 545 p. 


Manuscrito recibido el 30 de enero del 2006 
Aceptado el 18 de abril del 2008 

Este documento se debe citar como: 

Rodríguez-Laguna, R., J. Meza-Rangel, J. Vargas-Hernández y J. Jiménez-Pérez. 2009. Variación en la cobertura 
de suelo en un ensayo de procedencias de Pinus greggii Engelm. en el cerro El Potosí, Galeana, Nuevo 
León. Madera y Bosques 15(1):47-59. 



Página en blanco intencionalmente 
(This page left blank intentionally) 


Madera y Bosques 15(1), 2009:61-73 


61 


IARTICULO DE INVESTIGACION! 


Modelación del volumen fustal de Pinus 
durangensis en Guachochi, Chihuahua, México 

Modeling bole volume of Pinus durangensis in Guachochi, 

Chihuahua, México 

Marín Pompa-García 1 *, Ciro Hernández 2 , 

José Ángel Prieto-Ruiz 3 y Raymundo Dávalos Sotelo 4 

RESUMEN 

Pinus durangensis es la especie arbórea de mayor uso y distribución en la región de 
Guachochi, Chihuahua, México; sin embargo, el conocimiento para estimar con veracidad su volumen 
es limitado, por lo cual en este trabajo se planteó modelar el ahusamiento de dicha especie que posi- 
bilite generar una ecuación de volumen comercial compatible garantizando su aplicación operativa. 

Con ese propósito, se ajustó un modelo trigonométrico a 4 502 pares de datos de 663 árboles derri- 
bados en operaciones de aprovechamiento forestal. La prueba de bondad de ajuste consistió en la 
comparación de estadísticos como el coeficiente de determinación, la raíz del error medio cuadrático 
y un índice de precisión. El ajuste incluyó el uso de métodos estadísticos apropiados para detectar y 
corregir los problemas de autocorrelación y heterocedasticidad asociados al desarrollo de ecuaciones 
de ahusamiento. El modelo ajustado arrojó buenos resultados en su capacidad predictiva y mediante 
su integración matemática se obtuvo una ecuación de ahusamiento compatible con los volúmenes 
comercial y total. Se concluyó que las ecuaciones generadas tienen un sustento estadístico sólido, por 
lo cual son una herramienta fundamental para mejorar la precisión en los inventarios forestales, 
asegurando así su aplicación no sólo en la región de estudio, sino que también se recomienda su 
puesta en práctica en otras localidades, mediante las adaptaciones a las que haya lugar. 

PALABRAS CLAVE: 

Ahusamiento, Chihuahua, modelo trigonométrico, volumen. 


ABSTRACT 

Pinus durangensis is the most utilized and abundant tree species in the Guachochi región of the 
State of Chihuahua; however, there is a lack of appropriate knowledge to get accurate volume estima- 
tions. Therefore, the purpose of this study was to model the tapering of the bole to allow developing 
an equation to predict merchantable and total volume, which guarantees its operative application. For 
this purpose, a trigonometric model was fitted to 4 502 data pairs of 663 trees felled during harvesting 
operations. The goodness of fit test consisted of comparing the statistics, coefficient of determination, 
root mean square error, and a precisión índex. Model fitting ineludes the use of appropriate statistical 
procedures to account for the problems of autocorrelation and heterocedasticity associated with the 
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development of taper equations. The 
model gave best results in terms of 
predictive capacity, and via mathematic 
integration, an equation for tapering was 
obtained which is compatible with the 
merchantable and the total volume. It 
was concluded that the generated equa- 
tions have a solid statistical foundation, 
and therefore constitute an essential tool 
for improving accuracy in forestry inven- 
tories, assuring its usage not only in the 
study región but its application is recom- 
mended also in other sites, with the 
necessary adaptation of the model to the 
respective local conditions. 

PALABRAS CLAVE: 

Chihuahua, Taper, trigonometric model, 
volume. 


INTRODUCCIÓN 

En el estado de Chihuahua existen 
grandes superficies de bosques de pino 
cuyo aprovechamiento constituye la prin- 
cipal fuente de ingresos para los silvicul- 
tores. Particularmente, en el área de 
influencia de la organización Silvicultores 
Unidos de Guachochi, S.C., en el sur de 
Chihuahua, la importancia de Pinus duran- 
gensis es sobresaliente, pues se distribuye 
en 80% de la superficie arbolada y genera 
un volumen aprovechable anual del orden 
de los 190 000 m3 (SEMARNAT, 2005), el 
cual se industrializa para obtener 
productos de alto valor comercial, como 
madera contrachapada y aserrío, principal- 
mente. Sin embargo, es limitado el conoci- 
miento sobre la veracidad de las estima- 
ciones del volumen de esta especie, dado 
que en los programas de manejo de la 
región se utilizan ecuaciones desarro- 
lladas a nivel de género, obtenidas incluso 
para otras localidades y sin mayor 
sustento estadístico (SUG, 2007); específi- 
camente, según información de la direc- 
ción técnica, se ha venido utilizando el 
modelo de Biging (1984), que se ha 
empleado para la misma especie en 
regiones similares pero presenta desven- 
tajas en la estimación de diámetros para 


categorías mayores de 30 cm, sugirién- 
dose investigación adicional (Corral et al., 
1999). Por lo anterior, es imperativo contar 
con herramientas confiables para estimar 
volúmenes de madera comerciales y 
totales, a partir de ecuaciones que 
modelen el ahusamiento fustal; estas 
herramientas deben estar fundamentadas 
en procedimientos estadísticos rigurosos, 
para garantizar su aplicación operativa de 
acuerdo a las exigencias que marca la 
legislación forestal vigente. 

Dentro de la dasonomía existe un 
alto grado de desarrollo en la modelación. 
El empleo de ecuaciones de ahusamiento 
es un buen ejemplo de ello. A través del 
uso de metodologías adecuadas, particu- 
larmente de los procedimientos de regre- 
sión y de la modelación matemática, se 
han logrado considerables reducciones 
de costos en los inventarios, sin tener que 
disminuir la precisión de los estimadores 
en la búsqueda de optimizar la distribu- 
ción de productos forestales maderables 
previo a su utilización (Pompa y Solís, 
2008). 

Existen numerosos esfuerzos para 
modelar el ahusamiento de los fustes 
como apoyo para estimar su volumen 
comercial; en los últimos treinta años se 
han estudiado desde los modelos más 
simples de ahusamiento (Kozak et al., 
1969; Ormerod, 1973; Hilt, 1980), 
modelos polinomiales segmentados 
(Bruce et al., 1968; Max and Burkhart, 
1976; Cao etal., 1980) y modelos geomé- 
tricos y trigonométricos (Parresol y 
Tomas, 1996; Fang y Bailey, 1999; Zhang 
et al., 2002; Bi, 2000). Sin embargo, 
Newnham (1998) menciona que hay al 
menos dos razones para continuar con el 
estudio de este tema: 1) hasta hoy no 
existe una teoría que explique adecuada- 
mente la variación de la forma del fuste 
para todos los árboles y 2) las ecuaciones 
de ahusamiento pueden ser una herra- 
mienta eficiente para estimar tanto el 


Madera y Bosques 15(1), 2009:61-73 


63 


volumen comercial como el volumen total 
del fuste. Así, Clutter (1980) notó que en 
una ecuación dada de volumen comercial 
está intrínsecamente definida una función 
de ahusamiento, lo cual a su vez, 
mediante integración matemática, posibi- 
lita obtener una ecuación compatible de 
volumen total. 

Recientemente se ha demostrado la 
gran utilidad de usar modelos trigonomé- 
tricos (B¡, 2000; Bi y Long, 2001), ya que 
estos pueden ajustarse con gran precisión 
a datos procedentes de diferentes especies 
y formas, dado que el punto de inflexión 
varía de acuerdo al tamaño del árbol, situa- 
ción que se diferencia de los modelos 
típicos de forma o exponente variable. 

En general, los modelos de ahusa- 
miento deben cumplir con las siguientes 
características: 1) cuando la altura parcial 
es igual a la altura total, el diámetro 
predicho debe ser cero; 2) a 1 .3 m de altura 
debe predecir el diámetro normal y 
3) debajo de 1 .3 m de altura debe predecir 
diámetros mayores al diámetro normal. 


OBJETIVO 

El objetivo en este trabajo fue 
modelar el ahusamiento para Pinus 
durangensis, que permita, en conse- 
cuencia, generar ecuaciones de volumen 
comercial y total en los bosques de la 
organización Silvicultores Unidos de 
Guachochi, Chihuahua, México. 

MATERIALES Y MÉTODOS 
Descripción del área de estudio 

El área de estudio corresponde a la 
región administrativa de Silvicultores 
Unidos de Guachochi, S.C., cuya super- 
ficie es de 767 000 ha y se localiza en el 
suroeste del estado de Chihuahua, entre 
los paralelos 27°52’ y 27°30’ latitud norte 
y los meridianos 106°30’ y 107°30’ 
longitud oeste (Figura 1). Comprende una 
parte de los municipios de Urique, Bato- 
pilas, Balleza y Guachochi. 



Figura 1. Localización del área de estudio. 
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La Organización está integrada por 
30 ejidos y 1 000 pequeñas propiedades, 
con una población de 110 000 beneficia- 
rios, de los cuales 85% corresponde a la 
población Rarámuri (INEGI, 2004). Esta 
organización forestal se encuentra encla- 
vada en lo alto de la Sierra Madre Occi- 
dental, conocida como Sierra Tarahu- 
mara, con terrenos bastante acciden- 
tados que presentan elevaciones y depre- 
siones muy irregulares (Pompa, 1995). 

La vocación del uso del suelo es 
forestal; existe un potencial productivo de 
madera de 250 000 m 3 de pino y 80 000 m 3 
de encino por año (SEMARNAT, 2005). El 
bosque de esta región es de clima 
templado y semifrío con diferentes grados 
de humedad; predominan los géneros: 
Pinus, Quercus, Fraxinus y Juniperus, 
entre otros (SEMARNAT, 2005). 

Descripción de la población de estudio 

La población de estudio corres- 
pondió a las masas forestales arboladas 
donde vegeta Pinus durangensis dentro 
del área de influencia de Silvicultores 
Unidos de Guachochi, empresa de 
productores concesionaria de los servi- 
cios técnicos forestales, delimitándose 
mediante el uso de la cartografía dispo- 
nible en la dirección técnica sobre tipos 


de vegetación contenida en los 
programas de manejo de la zona. 

Descripción de la muestra 

Se tomó una muestra aleatoria de 
663 árboles para representar todas las 
calidades de estación del área de estudio, 
así como también la distribución de los 
árboles en términos de clases de 
diámetro y altura (Figura 2a). La toma de 
datos de campo se realizó mediante un 
muestreo destructivo aprovechando las 
áreas de corta del año 2003, para lo cual 
se midieron las siguientes variables: 
diámetro normal con corteza (D, en cm), 
altura total (H, en m), diámetro con 
corteza para cada sección ( d , en cm) y la 
altura de cada sección con respecto al 
suelo (h en cm). Se integraron 4 502 
pares de datos (d y h) distribuidos a lo 
largo del tronco. En la tabla 1 se muestra 
un resumen de las estadísticas descrip- 
tivas de los árboles utilizados. 

Donde: 

D = diámetro con corteza a la altura 
del pecho (1,3 m en cm). 

H = Altura total (m); h = altura de la 
sección del árbol a partir del suelo (en m). 

d = diámetro con corteza (en cm) a 
la altura h (en m). 


Tabla 1. Estadísticas descriptivas de los árboles utilizados. 


Variable 

Número de 
observaciones 

Media 

Desviación 

estándar 

Máximo 

Mínimo 

D 

663 

38 

10,7 

66 

10,2 

d 

4 502 

25,9 

16,2 

84 

0 

H 

663 

18,9 

5,8 

36,2 

4,2 

h 

4 502 

9,1 

7,75 

36,2 

0,1 
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Figura 2. Dispersión de observaciones diámetro normal - altura (a). Dispersión de 
alturas relativas respecto a diámetros relativos de los 663 árboles utilizados en el ajuste 

del modelo trigonométrico (b). 


La figura 2b representa las alturas 
relativas ( h/hi) frente a los diámetros rela- 
tivos ( d/D ) de cada sección medida del 
tronco. La amplitud de los datos refleja la 
magnitud de la forma de los árboles que 
conformaron la muestra utilizada. 


Modelo probado 

Para tratar de predecir el perfil de los 
árboles, con base en los datos prove- 
nientes de la muestra del área de estudio, 
se probó el modelo trigonométrico de Bi 
(2000). Su ecuación simplificada es 
dr = B K , donde dr es el diámetro relativo 
del fuste a una altura relativa dada, B es 
la base del modelo en función de la altura 


relativa y de las características dendro- 
métricas del árbol (diámetro normal y 
altura total), y K es una función que 
describe los cambios de la forma del fuste 
y sus características dendrométricas. La 
ecuación de base (B) ha sido construida a 
partir de ecuaciones trigonométricas de 
proporción del volumen según la geome- 
tría del fuste. Gracias a su flexibilidad, 
este modelo se ha utilizado para describir 
el perfil del fuste de Eucalyptus sp. en 
Australia (Bi, 2000), Pinus radiata en New 
South Wales (Bi y Long, 2001) y reciente- 
mente en Populus x euramericana en 
Navarra (Rodríguez y Molina, 2003), pero 
se desconocía su bondad de ajuste para 
Pinus durangensis. Su expresión mate- 
mática corresponde a: 
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d = D 
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Donde: 

d = diámetro con corteza que se 
alcanza a la altura h, en cm. 

h = altura desde la base del árbol 
hasta donde se alcanza el diámetro d, 
en m. 

D = diámetro normal con corteza, 
en cm. 

H = Altura total del árbol, en m. 

¡ = Coeficientes de regresión a deter- 
minar mediante el ajuste. 

Los cuatro primeros parámetros de 
la función K (fi 0 a P 3 ) describen los 
cambios de la forma del fuste a lo largo de 
éste, mientras que los tres últimos (p 4 a 
Pe) caracterizan la forma del fuste depen- 
diendo de las características dendromé- 
tricas del árbol evaluado. 

Este modelo cumple con la condi- 
ción de que d = 0 cuando H = h y pueden 
ser usados para predecir diámetros (d) a 
cualquier altura del trozo (/?), y volumen 
entre dos diámetros o alturas del tronco. 

Ajuste del modelo 

El ajuste del modelo se realizó por el 
método de mínimos cuadrados ordina- 
rios, empleando el programa estadístico 
SAS/STAT. (SAS Institute Inc., 2004). Al 
desarrollar procedimientos de regresión, 
usualmente se asume que el error es 
independiente y que se distribuye normal- 
mente con varianza y media igual a cero 


(Gujarati, 2004). Sin embargo, la hetero- 
cedasticidad y la autocorrelación son los 
inconvenientes más importantes en el 
caso de modelos de ahusamiento (Kozak, 
1997). Por tanto, en este trabajo se utili- 
zaron métodos generalizados de regre- 
sión lineal y no lineal para modelar la 
estructura del error. Los trabajos de 
Monserud (1984), Goelz y Burk (1992, 
1996, 1998) y Huang (1997) presentan 
aplicaciones prácticas de estos métodos. 
Los valores de la prueba de Park “P” 
(Park, 1966), y el estadístico de Durbin y 
Watson “DW” (Durbin y Watson, 1951), 
evidenciaron que la estimación final de 
los parámetros no estuvo afectada por 
problemas de heterocedasticidad y auto- 
correlación, respectivamente. 

Los criterios empleados para evaluar 
la bondad de ajuste del modelo se 
basaron en un análisis numérico y gráfico 
de los residuos, es decir, la diferencia 
entre el diámetro medido y el predicho. El 
análisis numérico consistió en la compa- 
ración de tres estadísticos utilizados con 
frecuencia en el ajuste de modelos: 1) el 
sesgo ( E ), que evalúa la desviación del 
modelo con respecto a los valores obser- 
vados; 2) la raíz del error medio cuadrá- 
tico ( REMC ), que analiza la precisión de 
las estimaciones, y 3) el coeficiente de 
determinación ajustado R 2 ac jj que repre- 
senta la parte de la varianza explicada 
por el modelo, teniendo en cuenta el 
número de parámetros de éste. Sus 
expresiones son: 
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Sesgo: 

Raíz del error medio cuadrático: 


Coeficiente de determinación ajustado: 


Donde y, , y, , y son el valor obser- 
vado, predicho y medio de la variable 
dependiente, respectivamente; n es el 
número total de observaciones usadas 
para ajustar el modelo; y p es el número 
de parámetros del modelo. 

Para el análisis gráfico se exami- 
naron diferentes representaciones de los 
datos experimentales y de los residuos 
después del ajuste de los modelos. Estos 
gráficos representan una herramienta 
muy importante en la selección de un 
modelo, ya que permiten detectar errores 
o comportamientos anormales (Rawlings, 
1988; Neter et ai, 1990). 

RESULTADOS 

En la tabla 2 se presentan los 
valores de los parámetros estimados y de 
los estadísticos del ajuste del modelo. 
Todas las estimaciones de los parámetros 
resultaron ser significativas al 95%. De 
acuerdo con los resultados mostrados en 
la tabla 2, el modelo explica de manera 
adecuada el perfil del árbol, ya que 
presenta valores bajos del REMC, expli- 
cando el 95,5% de la varianza del 
diámetro (d). 

Con relación a la heterocedasti- 
cidad, el gráfico de los residuos frente a 
los valores predichos por el modelo 
(Figura 3) y la prueba de Park (1966) 
evidenciaron que existe igualdad de 
varianza y que no hay relación estadísti- 
camente significativa de los residuales 


E = j?(.y ; -y.)ln 
REMC = 

R adj = 1 ' (” “ 1 ) ' 2 (y¡ ~ y¡ ) 2/(n - p )'2 (ji ~ y ‘ )2 
7 =1 7 = 1 


con las variables explicatorias. El estadís- 
tico de Durbin y Watson es próximo a 2 
(DW = 1 .3926) y está dentro de los inter- 
valos aceptados en las tablas respectivas 
para rechazar la existencia de autocorre- 
lación entre los residuos con un a de 0,05 
(Durbin y Watson, 1951). 

Aunque estos estadísticos y sus 
supuestos son buenos indicadores de la 
bondad global de ajuste de la función de 
ahusamiento, es conveniente una repre- 
sentación gráfica por secciones relativas 
de altura en el árbol. Según Castedo 
(2003), el modelado por posición es muy 
importante para decidir si en realidad la 
función de ahusamiento debe ser reco- 
mendada. Con esa finalidad, en la figura 
4 se ilustra el efecto que tiene el sesgo y 
el cuadrado medio del error a lo largo del 
fuste, para lo cual se gráfica la conducta 
del modelo. En los gráficos referidos, es 
evidente la bondad de ajuste del modelo, 
ya que tiene una mejor estimación del 
perfil en la parte basal donde se sitúa la 
mayor parte del volumen comercial; bajo 
este criterio de operatividad práctica se 
considera que este modelo es definitiva- 
mente un buen predictor de diámetros a 
diferentes alturas del fuste y por ende, un 
buen descriptor del ahusamiento de Pinos 
durangensis de la región de estudio. Su 
buen comportamiento estriba en su gran 
flexibilidad, ya que la altura normal rela- 
tiva (1.3/h) varía con la altura total del 
árbol y las características del exponente 
están relacionadas con la forma a lo largo 
del tronco y consideran las diferencias 
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Tabla 2. Comparación del ajuste de modelos de ahusamiento (con un nivel de 

significancia Pr > |t| < 0.0001). 



Figura 3. Gráfico de residuos frente a valores predichos del modelo de Bi (2000). 
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Figura 4. Comportamiento del sesgo y el REMC a lo largo del fuste para la función de 

ahusamiento de Bi (2000). 


entre árboles de distintos tamaños, lo 
cual es correspondiente a los resultados 
encontrados por varios autores para 
especies de pino (Castedo, 2003; Rodrí- 
guez y Molina, 2003). 

En la figura 5 se presenta la bondad 
de ajuste del modelo para los valores 
observados de c/a lo largo del fuste de un 
árbol seleccionado al azar. Se puede 
apreciar que el modelo simula satisfacto- 
riamente el perfil en la mayor parte del 
fuste, disminuyendo ligeramente su 
consistencia en la punta, donde se 
presenta un ahusamiento difícil de 
modelar debido a factores como el 
desarrollo de la copa, calidad del sitio y la 
densidad (Larson, 1963). 

Con base en los resultados discu- 
tidos, en la ecuación definida para la 


especie en estudio queda implícito el 
índice de esbeltez del fuste (H/D), 
mediante su inverso, es decir, D/H, como 
una medida de la proporcionalidad de 
esas dimensiones. Además, para cuando 
h=H, el diámetro estimado es cero, como 
lo indica la suma de sus coeficientes, 
cumpliendo satisfactoriamente con esta 
propiedad. 

DISCUSIÓN 

El modelo de Bi (2000) recomen- 
dado para Pinus durangensis en este 
estudio, tiene una importante ventaja 
comparativa con respecto al trabajo de 
Corral et al. (1999). En primer lugar, es 
consistente en sus estimaciones a lo 
largo del perfil del árbol, mientras que en 
la publicación referida se manifiesta la 
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perfil del árbol 


h (m) 



d observados 
Q Bi (2000) 


d(cm) 


Figura 5. Gráfico del ajuste al perfil de un árbol obtenido con el modelo de “Bi (2000)” 

propuesto para Pinus durangensis. 


presencia exponencial de errores para 
diámetros mayores a 30 cm con sus 
efectos multiplicativos en el cálculo 
erróneo del volumen, por lo cual los 
mismos autores recomiendan estudiar la 
constante de su modelo para mejorar el 
ajuste, el cual solamente contiene dos 
parámetros en su expresión, en contraste 
con el modelo aquí sugerido, cuya expre- 
sión es más compleja, aunque su integra- 
ción numérica no representa mayor difi- 
cultad para las modernas técnicas 
computacionales disponibles, como lo ha 
demostrado Castelo (2003) mediante el 
procedimiento QUAD del SAS/STAT® 
(SAS Institute Inc., 2004). 

Otra aportación importante del 
modelo es su compatibilidad con una 
ecuación de volumen (p. e., Smalian); es 
decir, cuando esta función es integrada 
genera estimaciones de volumen comer- 


cial y total idénticas a las obtenidas por 
una ecuación de volumen existente 
(Demaerschalk, 1972; Avery y Burkhart, 
1983). Sin embargo, debe considerarse 
que los volúmenes reales generalmente 
no se conocen, sino que suelen calcu- 
larse por las fórmulas de Smalian, Huber 
o Newton. 

CONCLUSIONES 

El modelo de Bi (2000) presentó 
buena capacidad predictiva para describir 
el perfil del fuste de Pinus durangensis en 
la parte suroeste del estado de 
Chihuahua. Este modelo puede ser usado 
confiablemente para predecir el diámetro 
a cualquier altura del fuste. Dicho modelo 
tiene la importante ventaja de ser una 
función de ahusamiento consistente en 
sus estimaciones a lo largo del perfil 
fustal; además, es compatible con las 
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ecuaciones de volumen existentes para el 
área de estudio. La ecuación generada 
en este trabajo constituye una herra- 
mienta valiosa para ser aplicada en los 
inventarios forestales en la región de 
estudio y tiene un sustento estadístico 
sólido, por lo que también se recomienda 
su aplicación en otras regiones similares, 
con las adaptaciones locales a las que 
haya lugar. 
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ARTICULO DE INVESTIGACION! 


Modelación matemática del secado 
convencional de madera 

Mathematical modeling of the conventional drying of wood 

Sadoth Sandoval-Torres* 


RESUMEN 

La operación de secado de madera es una etapa crucial dentro del procesamiento de productos 
de madera. Durante el secado se generan algunos defectos debido a la distribución variable del conte- 
nido de humedad y de la temperatura al interior del material. En el secado industrial de madera, los 
problemas de optimización de los procesos deben ser fundamentados en los diferentes fenómenos de 
transporte que controlan la transferencia de masa y calor en el sólido. Por sus propiedades físicas, las 
maderas duras adquieren gran importancia como materia prima, pues su dureza, color, estabilidad 
dimensional y alta resistencia a la degradación causada por el medio ambiente, les dan un valor 
comercial importante. En los últimos años el modelamiento de los procesos de secado ha sido una 
herramienta eficaz para la optimización de tales operaciones; estos modelos asocian balances de 
energía, de masa y de momento, que pueden ser expresados como ecuaciones diferenciales que son 
resueltas numéricamente para obtener soluciones específicas en cada problema. En este trabajo se 
explican los mecanismos de transferencia de masa y calor presentes en el secado convencional, se 
escriben las ecuaciones de transporte de masa y calor, y por último se muestran los perfiles de 
humedad y temperatura para el encino europeo Quercus pedonculae, sometido a secado conven- 
cional. 

PALABRAS CLAVE: 

Física de Secado, Transferencia de masa y calor, Madera latifoliada. 


ABSTRACT 

The kiln-drying of lumber is a very important stage into wood-based producís Processing. During 
drying process, some defects are produced due to the variable moisture and heat distribution inside of 
material. In order to optimize the kiln-drying of lumber is important to take into account the transport 
phenomena controlling the heat and mass transfer in the solid. Due to its physical properties, the hard- 
woods are a very important raw material, since its hardness, colour, dimensional stability and high 
resistance to weathering give them a very important commercial valué. At the last years, the drying 
process modelling has been an effective tool to optimize such operations; these models associates 
energy, mass and momentum balances, which can be expressed as differential equations that are 
numerically solved to obtain particularly solution for each problem. In this work, the heat and mass 
transfer mechanisms presents in kiln-drying of lumber are explained, heat and mass transport equa- 
tions are written, and finally moisture and temperature profiles for European Oakwood Quercus pedon- 
cula during kiln-drying are presented. 

KEY WORDS: 

Physics of Drying, Heat and Mass transfer, Hardwood 
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INTRODUCCIÓN 

El principal objetivo del secado 
industrial de madera es acelerar los 
mecanismos de transferencia de masa, 
optimizar el consumo energético y, princi- 
palmente obtener un producto con 
mejores atributos de calidad. El secado 
es una operación unitaria de uso 
frecuente en varias ramas de la industria. 
Durante una operación de secado, las 
maderas experimentan cambios físicos y 
químicos (Keey et al. 2000), los cuales 
afectarán su calidad al final del secado, 
es por eso que el estudio de la transfe- 
rencia de masa y calor en materiales 
biológicos adquiere gran importancia. De 
acuerdo a Martinovic et al. (2001), el 
secado es un paso muy importante dentro 
del proceso de manufactura de productos 
de madera, durante el cual se desarrolla 
una distribución de humedad y tempera- 
tura no uniforme en el material, que se 
refleja en muchas ocasiones en la gene- 
ración de defectos de secado. Los funda- 
mentos de los fenómenos de transporte 
que son aplicados al secado de madera 
pueden ser consultados en referencias 
como: Middleman (1997), Incropera et al 
(2006), Geankoplis (2003), Majumdar 
(2005) y Bird et al (2006). 

Los modelos matemáticos aplicados 
al secado de materiales están divididos 
principalmente en tres grupos: 1) modelos 
empíricos aplicados a procesos especí- 
ficos, 2) modelos difusivos de transfe- 
rencia de masa y calor que llevan a ecua- 
ciones simultáneas, y 3) modelos cono- 
cidos como comprensivos, los cuales 
involucran ecuaciones diferenciales de 
transporte de masa, calor y momentum 
con fundamentos termodinámicos; estos 
modelos son conocidos también como 
mecanísticos. En la actualidad el aporte 
científico del modelamiento matemático 
de los procesos de secado es de gran 
interés por su aplicación industrial, de ahí 
el interés de diversos centros de investi- 


gación por la simulación numérica de 
estas operaciones, la cual es una tarea 
de gran complejidad numérica y matemá- 
tica, como puede ser constatado en 
trabajos como el de A. Esfahanian (1999), 
Turner y Perré (2001) y Perré y Turner 
(2002) entre otros. 

OBJETIVOS 

1. Escribir las ecuaciones y explicar 
los mecanismos de transferencia de 
masa y calor presentes en el secado de la 
madera. 

2. Describir el secado de la madera 
como una operación unitaria fundamen- 
tada en los fenómenos de transporte apli- 
cados a un medio poroso. 

3. Presentar un modelo matemático 
mecanístico para el secado de madera y 
resolverlo numéricamente. 

Materiales porosos y materiales 
capilaro-porosos 

La distinción entre materiales porosos 
y materiales capilaro-porosos está basada 
en la presencia y tamaño de los poros. Los 
materiales porosos son definidos como 
aquellos que tienen un diámetro de poro 
mayor o igual a 10- 7 m, y los materiales 
capilaro-porosos tienen un diámetro menor 
a 10 7 m. El transporte de humedad en los 
materiales capilaro-porosos es un fenó- 
meno difícil de modelisar, pues además de 
la difusión molecular, el transporte de agua 
puede ser debido también a la difusión de 
vapor, difusión superficial, difusión 
knudsen, flujo capilar, flujo hidrodinámico 
entre otros mecanismos. 

En materiales no higroscópicos, los 
poros están llenos de agua si el material 
está saturado y con aire si ellos están 
completamente secos, estos materiales no 
se contraen. Por otra parte, en el caso de 
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materiales higroscópicos habrá un 
momento en el cual la presión de vapor de 
agua será importante, pues llega a ser 
función del contenido de humedad y de la 
temperatura del sistema. La madera perte- 
nece al grupo considerado como capilaro- 
poroso, el cual es higroscópico. En los 
materiales capilaro-porosos, los poros 
pueden estar comunicados, mientras que 
en los materiales porosos no necesaria- 
mente. Asimismo, la higroscopicidad en la 
madera, además de la porosidad, también 
es generada significativamente por su 
composición química, específicamente por 
la presencia de radicales OH. 

En el presente escrito, los valores y 
elementos de análisis están orientados al 
encino europeo la cual es una madera dura. 
De acuerdo a Williamson (2001), las 
maderas duras proceden de árboles de hoja 
ancha, pertenecientes a las angiospermas o 
conocidas como latifoliadas, y las maderas 
suaves a las procedentes de árboles con 
hojas aciculares, pertenecientes a las 
gimnospermas o confieras. Las maderas en 
forma general presentas dos secciones 
importantes, desde el punto de vista estruc- 
tural: una región conocida como albura y 
otra conocida como duramen. En la albura 
se llevará a cabo el transporte de nutrientes 
para el crecimiento del árbol, con el tiempo 
la porción central del tronco podrá colectar 
los excesos en nutrientes, los cuales meta- 
boliza para formar extractivos que otorgarán 
un color más intenso al duramen. La madera 
es un material poroso, biológico, anisotró- 
pico y heterogéneo (Truscott y Turner; 
2006), el cual posee una estructura fibrosa y 
celular muy compleja. 


réndales (MEF) de secado son basados en 
las leyes de transferencia de masa y calor, 
los cuales son resueltos con diferentes 
métodos numéricos, y llevan a soluciones 
particulares de cada problema. 

El flujo de contenido de humedad en 
el producto ocurre por difusión y es gober- 
nado por un coeficiente de difusión efec- 
tivo, siguiendo la ley de Fick, entonces el 
balance de masa en la tabla sera: 

D v-ÍVi o w)= ^ maclen; '^ 

^ humedad V v madera ' ) ^ ^ \ ) 


La ecuación anterior establece que 
la rapidez de cambio en el contenido de 
humedad W del material, es igual a la 
difusión de agua debido al gradiente de 
contenido de humedad interno de la tabla 


El producto se contrae durante el 
secado, sin embargo, los efectos de 
encogimiento, densidad y propiedades 
térmicas del producto son normalmente 
despreciados. 

En la actualidad, para fines prácticos 
y de modelación matemática, la cantidad 
de agua presente en la madera es clasifi- 
cada como agua libre, agua ligada y 
vapor de agua (Siau, 1984). El transporte 
de estos tres tipos de humedad son 
combinados para obtener un flux de agua 
simple j. Entonces la ecuación puede ser 
escrita como: 


Modelos difusivos 

En este trabajo sólo se describen los 
modelos difusivo y comprensivo para el 
secado de madera, pues los modelos 
empíricos no son de interés para este 
escrito. Los modelos de ecuaciones dife- 


J=D 


humedad 


V'(Vp, 


El flux másico de humedad j es 
proporcional a la difusividad del agua en 
el producto, D humedad 
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l 


P madera 


d Cp madera^ madera 


\ 


dt 


+ Cp 


dp 


mader 


madera J madera 


dt 


— = V • (k • VT madem ) - Cp agua j (y T madera ) (3) 


Ecuaciones para la transferencia de 
calor 


en el centro del material (X=0), por lo cual 
las condiciones límites en el centro son: 


La transferencia de calor al interior 
de una tabla de madera es dada por 
conducción, teniendo como parámetro 
motor el gradiente de temperatura que se 
desarrolla a lo largo del espesor de la 
tabla. La ecuación que describe el trans- 
porte de enegía en la tabla es escrita 
como se indica en la ecuación 3. 


En esa ecuación, el cambio de la 
energía interna en un volumen de control 
en un lapso de tiempo es representado 
por los dos términos a la izquierda de la 
ecuación. Este balance de energía esta- 
blece que la rapidez de cambio en la 
energía interna es igual a la diferencia 
entre la rapidez neta de calor transpor- 
tado hacia adentro del volumen de control 
debido a la conducción, 


V-{k 


vr 

madera madera 


J 


y la rapidez a la cual el calor es transpor- 
tado fuera del volumen de control debido 
al flujo másico, Cp agu J. 


Condiciones iniciales y límite 

Al inicio del proceso de secado la 
temperatura y el contenido de humedad 
del material son uniformes. Las condi- 
ciones iniciales son entonces: 


1,,.o - r„ 

(4) 

II 

O 

(5) 


dM 

dn 


n=0,t 


= 0 


dT 

dn 


n= 0/ 


= 0 


( 6 ) 

(7) 


En la superficie del material, la trans- 
ferencia de calor ocurre por convección, 
entonces podemos escribir: 

K(T, - T wt )- k ‘±L- j {*,„ +cp, m (T. ~ T. )]- 0 (8) 


El principio general de la convección 
en el secado es el siguiente: un flujo de 
calor superficial es debido a la diferencia 
que existe entre la temperatura del medio 
secante y la temperatura de bulbo 
húmedo (T x -T wb ). El flujo de calor que 
sale de la superficie de la tabla es cuanti- 
ficado por la conducción de calor en el 
producto y además por la energía interna 
del flux de humedad j que sale de la 
superficie. Cierta cantidad de calor será 
absorbida por el flujo de humedad. La 
ecuación 9 es aplicable cuando la 
humedad que está en la superficie de 
la tabla de madera no ha alcanzado su 
punto de evaporación y está todavía en 
forma líquida. Una vez que la humedad 
en la superficie de la madera es evapo- 
rada, el término de calor latente de 
evaporízación del agua hvap aglla , es elimi- 
nado y la condición límite en la superficie 
es escrita como: 


Debido a la simetría, no habrá un 
gradiente de humedad ni de temperatura 


K (r, -Tj-A-j-Cp^(T„-T,)- 0(9) 
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La transferencia de masa en la 
superficie ocurre por convección, la repre- 
senta la expresión dd 

dn 


que indica que el flujo de masa es normal 
a la superficie. La condición límite para la 
transferencia de masa convectiva en una 
superficie plana es: 


D, 


humedad 


dtPavJr) / \ ( 10 ) 

^ agua -= h p (w - w ) v ' 

m ' agua \ e oo / 


en donde h m es el coeficiente de transfe- 
rencia convectivo y W e es el contenido de 
humedad de equilibrio del producto bajo 
las condiciones específicas del aire 
(temperatura y humedad relativa, princi- 
palmente), asumiendo que el líquido y el 
vapor están en equilibrio en la superficie. 
Estas condiciones implican que la rapidez 
de eliminación de humedad debido a la 
convección es igual a la rapidez de flujo 
de humedad debido a la difusión. Esta 
condición límite se aplica al aire de 
secado, ya que el gradiente de contenido 
de humedad entre la superficie y el 
agente de secado es la fuerza motriz para 
la remoción de humedad del producto. 


transporte de masa y calor. El éxito de 
Whitaker (1977) fue el haber realizado un 
análisis riguroso de cómo transitar un 
sistema de un nivel microscópico a un 
nivel macroscópico con bases convin- 
centes. Las principales consideraciones 
en el desarrollo de Whitaker son un equi- 
librio termodinámico local, la aplicación 
de la ley de Darcy, la ley de Fick, un flujo 
de infiltración para el transporte de gas, 
un flujo capilar para el transporte de 
líquido, una estructura rígida del material 
(no hay encogimiento) y ausencia de 
agua ligada. La gran ventaja del modelo 
mecanístico es que la física del modelo 
es mejor entendida, las hipótesis son muy 
claras y los parámetros son bien defi- 
nidos. 

Se toma como hipótesis que existe 
un equilibrio termodinámico local, por lo 
cual las temperaturas promedio de las 
fases (sólida, líquida y gas) son iguales: 

_ S _ g _1 

T s =T g =T i = T (11) 


La presión de vapor en el medio es 
expresada como: 


_ g 

P v = aw Pv 

sat 


( 12 ) 


Modelo mecanístico para el transporte 

de masa en medio poroso p v = expí 46 784 - 6435 

1, T + 273.15 


( 13 ) 


De acuerdo a Ni et al. (1999), en 
teoría de secado, además de la proposi- 
ción fenomenológica de Luikov, Whitaker 
en 1977 usó un modelo mecanístico y 
desarrolló un grupo de ecuaciones para 
los medios porosos. Whitaker funda- 
mentó su modelo en las ecuaciones de 
conservación de masa y calor para cada 
fase presente en el material (sólido, 
líquido y vapor) y después de un 
promedio volumen de las fases, obtiene 
ecuaciones de continuidad y flux para el 


-3.8681n(r + 273.15/ 


Transferencia de masa 

A continuación se escriben las 
expresiones para la rapidez de transporte 
de las diferentes fases. 
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Para el líquido: 


p,v, 


k ■ k 

= Pl ^L-VP c -p, 
P¡ 


k k 

- =rl 

P¡ 


(vPg-p/g) 


\ / 


(14) 


Para el vapor: 


p v V v = p v —¡-^--S/Pg- p D-VC 
a 


(15) 


Se puede calcular un coeficiente de 
difusión de acuerdo a la proposición de 
Fuller, Schettler y Giddings (Poling et al., 
2000 ): 


io _ 3 r U5 


Dbin ~ 


M a g ua + M Aire 
M agua M Alre 






(20) 


Para el aire: 


Las ecuaciones de conservación de masa 


p n V a =p a ^2L-S7P s -p g D-VC 

P g ^ (16) 


El balance de masa se escribe para 
cada uno de los constituyentes en cada 
fase: 


La difusión-sorción toma la forma 
fenomenológica de un flux de agua 
ligada, la cual ha sido discutida abundan- 
temente en la literatura (Stamm, 1946; 
Stamm, 1960; Stanish et al., 1986): 


Para el aire seco y el vapor en la fase 
gaseosa se tienen: 


d p 
dt 


+ V- 


( - 8 - 

P,v 

\ 


\ 

= o 

/ 


( 21 ) 


J b = -p s D- VW b - p s D-W T ( 17 ) 

=b =bt 


d p 
dt 


+ V- 


( - 8 - \ 

P v V v 

V / 


K + K b 


( 22 ) 


Tranferencia de la fase gaseosa 


Para la fase líquida: 


A continuación se escriben las ecua- 
ciones para el transporte de la fase 
gaseosa (Shukla, 1990): 


d Pi 

dt 


+ V- 


/_ 1 \ 

p,y, 

V / 


= -K 


(23) 


Flujo viscoso de tipo Darcy: 


g 


k • k 



= = rg 


V. 


S7Ps 


Para el agua ligada: 

(18) 9 Pb + V • J =-K . 

dt b b 


(24) 


Difusión intermolecular de acuerdo a la 
ley de Fick 


J, ~-p g D‘ ff VC 

= V 

D eff = D v B 


(19) 


Entonces se puede obtener: 


dW 

dt 


■+ V- 


— (pi + p 8 K + Jb ) 
\P S ) 


= o 


(25) 


v 
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Realizando un balance enegético, de acuerdo con Whitaker (1977), se obtiene: 


dH 

~dT 


+ V' 


P, Vi Cp, +J b C Pl +Pa V a C Pa+ p* V v c Pv 


\ 


+ K Pv^ v+ K 


8 - \ 


PiVi + Pv v , 


-A •vr 


/ 


=0 


(26) 


Se establecen las condiciones límites para el modelo mecanístico: 


$ = 


í 


\ 


P¡Vi + p v V v+ J b 


•VA, 


(27) 


P¡V,+ p v V v + p b V b 


F. 


n = 


Ps 


(28) 


Pl y 1 Cp, + p b V b Cp, + p a Va Cp a + p v Vv Cp v 


T+ 


g _ 

I 

V 


K P v V v + h b p,V¡+p V 


8 - \ 


-a ■ vr 


/ 


\-n = Q + h° v+ h b+ C Pv T \F m (29) 


La presión capilar p c , como lo muestra la figura 1, es una función de la saturación 
del material (S) y de la temperatura, los efectos de la gravedad son ignorados. 



Figura 1. Presión capilar en encino europeo. 
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Modelamiento matemático del secado convencional de madera 


Se considera un contenido de 
humedad promedio de 0,9 kg agua/kg 
sólido seco. A continuación se presentan 
los parámetros que han sido determi- 
nados para madera de encino europeo 
(Hernández (1991): 


HR q = exp(M • b2 ]mw 


) 


(30) 


bl = 2.51x 10” 4 -T 2 -0.1811-7 + 35.719 


transformaciones de Boltzman (Diaw et al., 
2006) y de Laplace (Trofimov et al., 1992). 
Cuando la difusividad es una función de la 
humedad, lo cual puede ser en el caso del 
secado de la madera, las soluciones rigu- 
rosas no son viables. Para resolver una 
ecuación de difusión, la difusividad debe ser 
expresada en ciertas formas, por ejemplo 
como una expresión exponencial. Las apro- 
ximaciones adicionales son necesarias 
generalmente para simplificar las integrales. 


bl = -9.475 xlO” 4 -T +1.133 


Para la densidad básica de la madera: 


( 31 ) 


p(w)= p (0.626ff+l) 

i \ / i madera x / 


Presión capilar: 


P =1.5628x10 


H-Oexp(ifZ) 


(32) 


Coeficiente de difusión en el dominio 
higroscópico: 


A =■ 


1 


P madera 


-exp 


A 

T 


+a 2 \w + 



(33) 


A¡=10886.472 


A 2 =-3 1.232 


B ¡=-7600.705 


B 2 =6.697 


Soluciones analíticas y numéricas 

Las soluciones analíticas pueden ser 
rigurosas o aproximadas. Las soluciones 
analíticas rigurosas son solamente para 
problemas limitados como la ecuación de 
difusión simple con coeficientes constantes 
y geometrías regulares; algunos utilizan las 


Trabajos recientes como el de 
Babaszak y Kowalski (1999), presentan solu- 
ciones para la deformación de materiales 
durante el secado haciendo uso de las trans- 
formadas de Laplace y la separación de 
variables. Las soluciones analíticas poseen 
fuertes requerimientos para la formulación 
matemática, incluyendo ecuaciones diferen- 
ciales, coeficientes, términos fuente, condi- 
ciones límite y condiciones iniciales. 

En el caso de soluciones numéricas, 
es común en problemas de transporte de 
masa y calor, hacer uso de los métodos 
de elemento finito, volumen finito y dife- 
rencias finitas. El problema de modelar 
mecanismos en bioproductos, como el 
secado de madera, es que la solución 
depende de las propiedades, pero 
también de un fuerte acoplamiento de los 
procesos, dichos acoplamientos son en 
general no lineales. Ecuaciones 
acopladas de transferencia de masa y 
calor pueden ser resueltas utilizando 
esquemas iterativos o soluciones directas, 
dependiendo de la memoria disponible en 
el computador. Con respecto a la repre- 
sentación gráfica de resultados numé- 
ricos, en la actualidad es común el uso de 
programas comerciales que resuelven 
ecuaciones diferenciales, o que calculan 
flujos en medios porosos, por ejemplo 
FLUEN Inc., AEA-CFDS Inc., MATHCAD 
y MATLAB. Kocaefe et al. (2006), utilizan 
FEMLAB para simular el secado de 
madera aspen a alta temperatura. 
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RESULTADOS Y CONCLUSIONES 

En esta última parte se presentan los 
perfiles de contenido de humedad y de 
temperaturas para una tabla de madera de 
una pulgada de espesor, utilizando el 
modelo mecanístico utilizando mathcad. 
La madera es sometida a secado convec- 
tivo a una temperatura de bulbo seco de 
80°C, una humedad relativa de 30% y 
velocidad de flujo de aire de 3 m/s. Con 
respecto a la parte de resolución del 
problema numérico, se utiliza el método 
implícito en diferencias finitas para discre- 
tizar el sistema de ecuaciones diferen- 
ciales. Para resolver eficientemente el 
sistema de ecuaciones resultante, se 
aplica la factorización LU y el método 
Crank Nicholson para evitar problemas de 
inestabilidad (Baehr y Stephan, 1998). La 
finalidad de esta sección es mostrar al 
lector, de forma gráfica, que en la actua- 
lidad se pueden estimar los perfiles de 
humedad y temperatura en sólidos, consi- 
derando un modelo de secado completo, 
tomando en cuenta los fenómenos de 
transporte mencionados anteriormente. 
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En las figuras 2, 3, 4, y 5 se observan 
los perfiles de distribución de humedad en 
una tabla de madera para diferentes 
tiempos de secado. Como es de esperarse, 
los perfiles de humedad en general serán 
parabólicos, el transporte es lento, pues por 
una parte, la permeabilidad en maderas 
duras como el encino retrasan el transporte 
de humedad en el material, y por otra, la 
remoción de humedad en el dominio higros- 
cópico es lento, es decir, por abajo del punto 
de saturación. 

La figura 6 muestra los perfiles de 
temperatura en la tabla; en dicha figura se 
puede observar que la temperatura de la 
superficie es menor que la temperatura del 
aire (80°C) durante casi todo el proceso 
debido a la evaporación existente en la 
superficie del material; la evaporación 
provoca un enfriamiento de la superficie. Sin 
embargo, al final del proceso, es decir, 
cuando la tasa de evaporación es casi cero, 
la temperatura de la pieza será casi igual a 
la del aire, pues la cantidad de calor cedido 
al agua es menor, pues es menor la cantidad 
de masa de agua evaporada. 



Figura 2. Distribución de humedad en la tabla (2 horas). 
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Figura 3. Distribución de humedad en la tabla (10 horas). 
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Figura 4. Distribución de humedad en la tabla (48 horas). 
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Espesor de 


Figura 5. Distribución de humedad en la tabla (100 horas). 



Figura 6. Perfil de temperaturas en la tabla de madera. 
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Modelamiento matemático del secado convencional de madera 


La importancia práctica de predecir la 
distribución de humedad y temperatura en 
la madera, es conocer la manera en que se 
desarrollan los gradientes al interior del 
material, además de predecir en qué 
momento la diferencia de contenido de 
humedad entre el centro y la superficie es 
del 3%, pues de esta forma se reducen los 
esfuerzos de tensión que se generan, mini- 
mizando los riesgos de colapso, de apana- 
lamiento (grietas internas) y de las grietas y 
rajaduras superficiales. Al término de un 
proceso de secado, la presencia de 
esfuerzos debidos a los grandes gradientes 
de humedad al interior del material o en la 
superficie, pueden generar defectos de 
secado como son los colapsos, las raja- 
duras y torceduras, que disminuyen el valor 
comercial de las maderas. 

Con respecto a la temperatura, la 
importancia de predecir los perfiles de 
transferencia de calor, es que en general, 
en el secado de materiales biológicos, se 

Notación 


debe respetar una temperatura máxima 
en todo el material, de otra forma se corre 
el riesgo de inducir cambios en su estruc- 
tura física o bien degradar determinada 
composición química que es deseable en 
el producto final; de esta forma se puede 
tener una idea más clara del momento en 
que el centro del material ha alcanzado 
los requerimientos energéticos impuestos. 
En este caso, el perfil de temperatura 
presentado en la figura 6, muestra cómo 
se desarrollan los gradientes térmicos al 
interior de la madera. Esto es útil para 
optimizar los consumos energéticos del 
proceso; de esta forma se puede saber 
en qué etapa del secado se puede sumi- 
nistrar menor temperatura, o hacer 
cambio de condiciones del aire, ya que en 
la práctica del secado las operaciones 
industriales son conducidas por 
programas o secuelas de secado, que 
son una serie de etapas de proceso, y en 
cada etapa las condiciones de secado 
son diferentes. 


D=H Coeficiente de difusión efectivo (m 2 /s) 

D v Coeficiente de difusión de vapor en el aire (m 2 /s) 

k Permeabilidad intrínseca 

00 

Pg Densidad del vapor en la fase gaseosa (kg/m 3 ) 

Densidad del aire en la fase gaseosa (kg/m 3 ) 

D Difusividad equivalente (m 2 /s) 

k Permeabilidad relativa al líquido 

k g Permeabilidad relativa al gas 

P 9 V Presión promedio del vapor en el gas (Pa) 

Df f Coeficiente de Difusión efectivo del vapor (m 2 /s) 

e Porosidad 

p a ¡ r Viscosidad del agua líquida (Pa s) 

p g Viscosidad de la fase gaseosa (Pa s) 

v g Viscosidad cinemática del gas 

Pi Densidad del líquido (kg/m 3 ) 

Densidad básica de la madera (kg/m 3 ) 


Pmadera 
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v 

dw 

B 

c 

CPagua 
CP madera 

D b 

Dbin 

D h 

D humedad 
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hr 
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'm 


eq 


HR, 
hVggua 
J 

íb 

K 

K b 

km adera 

Magua 
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n 

p 

Pe 
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R 

S 

t 
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T s 

Twb 

W 

Woo 
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Volúmenes de difusión 
Actividad del agua 
Factor freno 

Concentración de vapor en el aire (kg vapor agua / kg aire seco) 
Calor específico del agua (J/kg °C) 

Calor específico de la madera (J/kg °C) 

Coeficiente de difusión del agua ligada (m 2 /s) 

Coeficiente de difusión binario (m 2 /s) 

Coeficiente de difusión en el dominio higroscópico 
Coeficiente de difusión del agua (m 2 /s) 

Aceleración de la gravedad (m/s2) 

Coeficiente de transferencia de calor (watt/m 2 °C) 

Coeficiente de transferencia de masa (m/s) 

Humedad relativa de equilibrio (%) 

Calor latente de evaporación del agua (J/kg) 

Flujo másico de humedad (kg/m 2 sec) 

Flujo másico de agua ligada (kg/m 2 sec) 

Tasa de evaporación del agua libre (kg agua libre/sec) 

Tasa de evaporación del agua ligada (kg agua ligada/sec) 
Conductividad térmica de la madera (watt/m °C) 

Peso molecular del agua (mol/g mol) 

Peso molecular del aire (mol/g mol) 

Normal 

Presión del sistema (Pa) 

Presión capilar (Pa) 

Presión de vapor saturado (Pa) 

Constante de los gases ideales 
Saturación líquida 
Tiempo (sec) 

Temperatura del aire (°C) 

Temperatura promedio de la madera (°C) 

Temperatura superficial del sólido (°C) 

Temperatura de bulbo húmedo (°C) 

Contenido de humedad de la madera (kg agua/kg madera seca) 
Contenido de humedad crítico (kg agua/kg madera seca) 
Contenido de humedad de equilibrio (kg agua/kg madera seca) 
Punto de saturación de la fibra (kg agua/kg madera seca) 
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Consideraciones geométricas 
en la determinación de las propiedades 
en flexión estática de bambú 

Geometric considerations for determining static bending 

properties of the bamboo 

Víctor Rubén Ordóñez Candelaria 1 
e Ignacio Salomón Quintana 2 

RESUMEN 

Se presenta un método alternativo para calcular las propiedades geométricas de las secciones 
de los culmos de bambú que se utilizan en ensayos de flexión para obtener las propiedades mecá- 
nicas bajo esta solicitación. El método de prueba normalizado recomienda se calcule el momento de 
inercia con el promedio de los diámetros y espesores de pared medidos en una sección, pero debido 
a las irregularidades naturales del culmo se han obtenido errores, con las secciones así medidas, de 
hasta 10%. Se propone un procedimiento utilizando imágenes digitalizadas de las secciones de los 
culmos que se procesan con un programa de CAD (Computer Aided Design), efectuando una integra- 
ción numérica, con este método el error se puede reducir a un valor aproximado de 0,011%. También 
se analiza otra fuente de error al calcular el MOE, cuando se usa corno valor máximo de la deforma- 
ción la que ocurre al centro del claro, pero por la variación de la sección la máxima deformación no 
se presenta en ese punto. 

PALABRAS CLAVE: 

Culmos, flexión, módulo de elasticidad, módulo de ruptura, momento de inercia. 


ABSTRACT 

ln this paper an alternative method to evalúate the geometric properties of the cross section of 
bamboo culms is presented, which are necessary to obtain the mechanical properties of this material 
subjected to bending tests. The standard test methods recommend calculating the second moment of 
area for a cross section of bamboo culms with the average diameter of the specimen, measuring two 
cross sections near of load polnts. But the section irregularities of culms yield errors of up to 10%. In 
this paper a procedure is presented to make a numerical integration using Computer Aided Design 
(CAD) software. Basically the procedure proposes to obtain a scanned image of the cross sectional 
area and then to digitalize it for use with the CAD software and carry out the numerical integration to 
get the moment of inertia. Using this method the errors are reduced to around 0,011%. Another 
possible source of error is analyzed when the MOE is calculated using the deflection valúes measured 
at mid-span, instead of the máximum valué which occurs outside of this point. 

KEY WORDS: 

Culms, modulus of elasticity, modulus of rupture, second moment of area, static bending. 
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Consideracionens geométricas en i a determinación de las propiedades en flexión... 


INTRODUCCIÓN 

La determinación de las propie- 
dades mecánicas de materiales que 
tienen aplicaciones estructurales requiere 
de cierto grado de confiabilidad para 
obtener los valores de esfuerzos básicos 
para el diseño de estructuras. Las 
pruebas para obtener información de las 
propiedades mecánicas de los mate- 
riales, en general, se realizan con espe- 
címenes que se elaboran con formas 
regulares y tolerancias muy estrictas en 
sus dimensiones y formas. Sin embargo, 
en el caso de algunos materiales fores- 
tales, es conveniente hacer pruebas 
mecánicas de resistencia sin modificar 
sus formas naturales, y en este caso no 
se tiene control sobre su geometría y sus 
dimensiones. En el caso de los tallos de 
bambú, la determinación de las propie- 
dades mecánicas en flexión del culmo 
requiere trabajar con secciones cercanas 
a la circular, aunque sus irregularidades 
pueden llegar a ser muy grandes. La 
norma ISO 22157 (2004), que se aplica 
para realizar pruebas mecánicas con 
culmos de bambú, establece que se 
calculen sus propiedades geométricas 
con el promedio de diámetros y espe- 
sores de la pared de una sección. En el 
caso de flexión se miden las secciones 
en los extremos del espécimen, para 
control de la prueba, al final de la misma 
se miden las secciones cercanas a los 
puntos de carga, y con el promedio de 
estas mediciones se calculan el módulo 
de ruptura (MOR) y el módulo de elasti- 
cidad (MOE). Así, estas variaciones 
geométricas son compensadas y el 
elemento es tratado (en los cálculos) 
como una viga de sección transversal 
perfectamente circular y constante a lo 
largo del tallo. Este método da una buena 
aproximación al comportamiento real de 
un elemento. No obstante, existen ya 
herramientas computacionales con las 
que es posible describir y manejar de 


mejor manera estas irregularidades 
geométricas, permitiendo que la determi- 
nación de las propiedades mecánicas 
sea más precisa. 

Durante la aplicación de esta norma, 
en un lote de pruebas para determinar las 
propiedades en flexión de los tallos de 
bambú, se detectó que los valores de 
MOR y MOE son considerablemente dife- 
rentes al obtenerlos con las recomenda- 
ciones de la norma y calcularlos con las 
propiedades geométricas reales de la 
sección, es decir, tomando en cuenta la 
irregularidad de la sección considerada 
circular. Para estimar la magnitud del 
error utilizando el método que describe la 
norma ISO 22157 contra los valores obte- 
nidos al considerar la sección real, se 
analizaron los valores de una muestra 
pequeña de un lote de pruebas en flexión, 
aplicando ambos procedimientos. Lo más 
importante de este trabajo es la presenta- 
ción de un procedimiento sencillo y de 
mayor exactitud para calcular el momento 
de inercia de una sección irregular, utili- 
zando imágenes digitalizadas que se 
procesan en un programa de CAD, y una 
estimación del error que se obtiene con 
este procedimiento. 


OBJETIVOS 

Desarrollar un procedimiento para 
calcular el momento de inercia de una 
sección irregular, utilizando imágenes 
digitalizadas y un programa de CAD, 
evaluando su exactitud. 

Estimar la magnitud del error al 
calcular el momento de inercia de las 
secciones transversales de los culmos, 
cuando se considera para su determina- 
ción que, la sección es circular y cons- 
tante a lo largo del espécimen, con 
respecto al momento de inercia obtenido 
mediante imágenes digitalizadas. 
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MATERIALES Y MÉTODOS 

Para este trabajo se utilizó una parte 
de los datos de un lote de ensayos de 
flexión con culmos de bambú con aproxi- 
madamente 110 especímenes. Las 
pruebas se realizaron conforme a lo que 
establece la norma ISO 22157-1 (2004). 
Se seleccionaron culmos sin defectos y 
con una longitud mínima de 30 veces su 
diámetro. La prueba se efectuó tal como 
lo establece la norma mencionada, apli- 
cando carga en cuatro puntos (figura 1), 
utilizando los aditamentos especiales en 
los apoyos y en los puntos de aplicación 
de la carga. La aplicación de la carga se 
realizó a una velocidad uniforme, determi- 
nando la carga máxima que resistió la 
pieza hasta la falla; se obtuvo una gráfica 
de carga-deflexión durante la prueba. Al 
finalizar la prueba se midieron las 
secciones transversales en tres puntos a 
lo largo del culmo, dos próximos a los 
apoyos y uno en el centro del claro de 
carga. Se midieron los diámetros y 
grosores de pared en 4 puntos de la 
sección transversal (figura 2a derecha). 
El valor promedio de estas mediciones se 
usó para calcular el momento de inercia 
(I) de sección transversal como si fuera 
una sección circular perfecta. 

Para calcular el módulo de ruptura 
(MOR) y el módulo de elasticidad (MOE) 


de la pieza ensayada, se usaron las 
siguientes ecuaciones: 

MOR = -^¿— 0) 

61 2 


Donde: 


MOR = Módulo de ruptura, MPa 
Pmax = Carga máxima que resiste el 
culmo, N 

L = Claro de carga o distancia entre 
apoyos, mm 

D = Diámetro externo del culmo en 
el punto central del claro de carga, mm 
I = Momento de inercia de la sección 
transversal del culmo en la parte central 
del claro de carga, mm 4 


MOE 


23PL 3 

129616 


( 2 ) 


Donde: 

MOE = Módulo de elasticidad, MPa 
P = Carga aplicada al espécimen 
en un punto del comportamiento lineal del 
material, N 

I = Momento de inercia promedio de 
la sección transversal del culmo, mm 4 
A = Deformación del espécimen al 
centro del claro, correspondiente a la 
carga P, mm 


carga 




viga repartidora de carga U 







-aplicador 
apoyo 


L 

3 


L 

3 



Figura 1. Esquema del ensayo de flexión estática con culmos de bambú. 
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Figura 2. a) Imagen digitalizada de la sección transversal de un culmo (izq), mediciones 
de los diámetros y grosores de la pared (der). b) propiedades geométricas de la misma 

sección obtenidas con el programa de CAD. 


La medición de la carga (P), la defor- 
mación del espécimen (A) y el claro de 
carga (L), en la prueba de flexión estática, 
se realiza con una exactitud aceptable. 
Sin embargo, la obtención de las propie- 
dades geométricas, y en particular el 
momento de inercia de una sección que 
se aleja considerablemente de un círculo 
no se realiza con la misma calidad que se 
tiene en los otros parámetros de las ecua- 
ciones 1 y 2. Una mejor manera para 
calcular el momento de inercia es 
mediante una integración numérica, sin 
embargo, en la práctica resultaría muy 
laborioso efectuar estos cálculos para 
cada sección. 


Se propone un método que consiste 
en digitalizar la sección del tallo más 
cercana a la falla (el momento de inercia 
de esta sección se usa para calcular el 
MOR) que no haya sido alterada durante 
la prueba. La imagen digitalizada con un 
escáner se exporta al programa Autocad; 
(Autodesk); sobre la imagen (formato jpg) 
se dibuja un esquema de la forma trans- 
versal del tallo siguiendo el contorno de la 
misma, utilizando las herramientas del 
programa de Autocad, se debe asegurar 
que el plano de aplicación de las cargas 
coincida con el eje vertical de la imagen 
en la figura digitalizada. Para calcular el 
Módulo de Elasticidad (MOE) se digitali- 
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zaron además la sección central más 
cercana a la falla y las secciones 
cercanas a los apoyos del espécimen 
ensayado. 

Con las imágenes digitalizadas se 
dibuja con el programa de Autocad su 
contorno para obtener la sección del 
espécimen. Una vez que se tiene la 
sección se pueden obtener las propie- 
dades geométricas de la sección dibujada 
utilizando el comando MASSPROP del 
Autocad, la figura 2a muestra la imagen 
real de la sección y en la figura 2b se 
observa la sección dibujada con el 
Autocad con su salida típica de datos del 
programa. Tomando una muestra de 13 
secciones de culmos se obtuvieron las 
diferencias entre el momento de inercia 
obtenido con la sección digitalizada y el 
calculado con el promedio de los diáme- 
tros y grosores, en la figura 3 se muestran 
los resultados obtenidos, se observan 
diferencias hasta del 10%, mismo valor 
que repercute directamente en la estima- 
ción de las propiedades mecánicas del 
bambú. 

Con el fin de evaluar la exactitud del 
procedimiento que utiliza la sección digi- 
talizada se realizó una estimación de los 
errores que se obtienen entre los 
momentos de inercia obtenidos con la 
sección digitalizada de un disco metálico 
fabricado en un torno y los valores que se 
obtienen midiendo el mismo disco con un 
calibrador (‘pie de rey’) digital que tiene 
una resolución de 0,01 mm, un error 
máximo de 0,02 mm y una incertidumbre 
de ± 9 pm. La forma en que se fabricó el 
disco metálico en un torno de metal- 
mecánica permite tener un círculo más 
cercano al ideal. En este ejercicio se 
midió ocho veces el diámetro del disco de 
acero en diferentes direcciones y se 
calcularon los momentos de inercia 
correspondientes. De la misma sección 
se obtuvieron ocho digitalizaciones con 
las cuales se obtuvieron los valores del 


momento de inercia utilizando las herra- 
mientas del programa de “AutoCad”, al 
comparar ambos valores se obtuvo una 
diferencia promedio de -0,011% que 
representa un valor muy inferior a los que 
han sido obtenidos para las secciones de 
bambú. Estos resultados se obtuvieron 
siguiendo el mismo procedimiento utili- 
zado con las secciones de los culmos. 

Otra fuente de error se tiene al medir 
las deformaciones de la viga al centro del 
claro, pero en realidad suceden en otro 
punto cerca del centro, para evaluar el 
error se realizó el siguiente ejercicio. Las 
deflexiones de una viga se obtienen 
resolviendo la ecuación diferencial de la 
curva de deflexión (3); en este trabajo se 
presentan las ecuaciones aplicables al 
tercio medio de la viga (Gere, 2002): 


El 


X 



(3) 


De donde: 


dv 2 

ck 2 


=MEI X 


(4) 


integrando una vez: 

ecuación de la pendiente de la viga 

dv M | c ( 5 ) 
dx El x 


integrando la ecuación de la pendiente: 
ecuación de la deflexión de la viga 


dv M ^ 

= ^ l 

dx El x 


( 6 ) 


Donde: 

M = PL = momento flexionante en el 
6 tercio medio de la viga 
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Figura 3. Diferencia entre el valor del momento de inercia con la sección digitalizada y 
el calculado con los diámetros y grosores promedio. 


E = Módulo de elasticidad en flexión 
del bambú = MOE 

l x = Momento de inercia de la 
sección del bambú, considerando varia- 
ción lineal del diámetro se calcula con la 
siguiente expresión 

P = Carga aplicada en la viga 
L = Claro de carga 



¿ 

D ü-( D „- D , 2 )f) 


( 7 ) 


D e i, D e 2 = diámetro exterior del 
culmo en los apoyos 1 y 2 

D¡i, D ¡2 = diámetro interior del culmo 
en los apoyos 1 y 2 


Siempre y cuando se conozcan las 
funciones del momento de inercia de la 
sección y de la variación del momento 
flexionante a lo largo de la viga, se puede 
determinar el punto donde se presenta la 
máxima deformación de la viga; la función 
que se debe integrar, suponiendo que la 
variación del diámetro y del espesor de la 
pared es lineal entre los extremos del 
culmo, la mostrada en la ecuación 8. 

Integrar la ecuación 5 representa un 
procedimiento complejo; para resolver 
este problema se puede acudir a la teoría 
de la resistencia de materiales que se 
tiene en la literatura correspondiente 
(Gere, 2002). Una manera más práctica 
de obtener las deformaciones de una viga 
de sección variable se logra aplicando el 
teorema de Área-Momento, el cual se 
aplica para determinar el punto exacto de 
la máxima deformación de una viga con 
sección transversal variable. 
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M 



( 8 ) 


En la figura 4 se presenta un 
esquema de cómo se obtendrían las 
deformaciones de la viga; se hizo el ejer- 
cicio con los valores de un espécimen 
ensayado en laboratorio con una varia- 
ción alta de su sección transversal. 


CARGAS Y DIAGRAMA DE MOMENTOS. PRUEBAS DE FLEXIÓN 




CALCULO DE LA DESVIACION TANGENCI ALtc, Y t» 



Teorema 1. El ángulo iB/A entre las 
tangentes a la curva de deflexión en dos 
puntos Ay B es igual al área del diagrama 
de M/EI entre estos dos puntos. 

Teorema 2. La desviación tangencial 
t B /A del punto B desde la tangente en el 
punto A es igual al primer momento está- 
tico del área del diagrama M/EI entre A y 
B, evaluado respecto a B. 

Haciendo el ejercicio se obtiene la 
ubicación del punto donde se presenta la 
máxima deformación de la viga que se 
localiza a 3,14 cm del centro. El error 
entre la deformación que se presenta a 
esta distancia contra la deformación 
medida en el centro del claro, tal como 
se realizó en los ensayos, es despre- 
ciable. 


RESULTADOS 

De la revisión de los experimentos 
realizados se observa que, medir las 
propiedades geométricas de las secciones 
con un programa de dibujo sobre una 
figura digitalizada arroja resultados con 
más exactitud que la propuesta de norma 
de tomar los promedios de las dimen- 
siones (Tabla 1). 

Se presentan los valores de las 
propiedades geométricas de una sección 
circular (rondana de acero), obtenidas 
tanto de las secciones digitalizadas como 
de las obtenidas tomando mediciones con 
vernier (Tabla 2). 


Figura 4. Cálculo de las desviaciones 
tangenciales tBA y tCA. 
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Tabla 1 . Valores de mediciones y momentos de inercia calculado y obtenido de 

secciones de culmo digitalizadas. 

sección D prom 
mm 

tprom 

mm 

Icalc 

mm 

Iscan 

mm 

Error 

relativo 

1 

71,79 

16,94 

1 202 457 

1 266 100 

-5,1% 

2 

111,17 

17,69 

5 877 930 

6 145 000 

-4,4% 

3 

80,68 

26,08 

2 047 380 

2 108 200 

-2,9% 

4 

55,83 

13,56 

443 564 

466 900 

-5,0% 

5 

54,86 

11,22 

390 396 

379 900 

2,7% 

6 

58,34 

12,01 

500 536 

540 500 

-7,4% 

7 

54,78 

11,83 

395 997 

401 300 

-1,4% 

8 

81,68 

11,62 

1 612 360 

1 465 300 

10,1% 

9 

77,58 

27,84 

1 766 861 

1 696 700 

4,1% 

10 

82,42 

11,77 

1 675 182 

1 717 900 

-2,5% 

11 

84,08 

14,53 

2 003 410 

2 020 400 

-0,8% 

12 

77,64 

23,74 

1 743 044 

1 828 600 

-4,7% 

13 

84,34 

27,66 

2 448 915 

2 509 400 

-2,4% 
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Tabla 2. Valores de momentos de inercia de una pieza metálica (rondana) de sección 
circular, obtenidos usando mediciones con vernier y calculados con la sección 

digitalizada. 


repeticiones 

limag 

Imed 

Error 

error relativo 

1 

1 502 564,68 

1 502 431,20 

133,48 

0,009 % 

2 

1 503 436,50 

1 502 431,20 

1 005,30 

0,067 % 

3 

1 501 406,60 

1 504 047,81 

-2 641,21 

-0,176 % 

4 

1 501 962,52 

1 502 431,20 

- 468,68 

-0,031 % 

5 

1 503 334,79 

1 501 623,39 

1 711,40 

0,114% 

6 

1 500 622,42 

1 500 815,90 

- 193,48 

-0,013% 

7 

1 501 803,15 

1 502 431,20 

- 628,06 

-0,042 % 

8 

1 502 161,52 

1 502 431,20 

- 269,68 

-0,018 % 

promedio 

1 502 161,52 

1 502 330,39 

- 168,87 

-0,011 % 

desviación estándar 

885,40 

851,02 



cv 

0,059 % 

0,057 % 




DISCUSIÓN 

Es indudable que la sección de los 
culmos del bambú tiene una gran varia- 
ción en sus dimensiones a lo largo del 
mismo, así como que la forma del mismo 
dista mucho de ser un círculo perfecto; 
esto último puede causar grandes errores 
en el cálculo de las propiedades mecá- 
nicas del bambú, utilizando los valores de 
las propiedades geométricas como lo 
establece la norma ISO-22157-1. Inde- 
pendientemente de que en la práctica el 
diseño no se realice con tanto refina- 
miento al momento de dimensionar 
estructuras con bambú, en este caso ya 
se tendrán otras consideraciones para el 
manejo de las secciones de los culmos. 


Considerando los resultados de este 
trabajo, la mejor forma de calcular las 
propiedades geométricas de secciones 
irregulares es determinarlas utilizando las 
herramientas que tiene el programa 
Autocad. Este procedimiento requiere de 
un poco más de trabajo, pero se logra una 
mejor estimación de la resistencia del 
bambú. Los errores obtenidos mediante 
este procedimiento tienen valor promedio 
de 0,011%, aunque éste puede mejorarse 
dependiendo de la exactitud con que se 
digitalice la imagen. Estos errores no son 
comparables con los obtenidos al usar el 
otro procedimiento. 
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CONCLUSIONES 

Los resultados de este trabajo 
permiten concluir que el uso de las 
secciones digitalizadas de los culmos de 
bambú, para determinar sus propiedades 
geométricas, genera mejores resultados 
de las pruebas que el procedimiento que 
propone la norma ISO 22157-1 (2004). El 
uso actualmente tan común de las herra- 
mientas de cómputo facilita este procedi- 
miento, aunque se requiera un poco más 
de trabajo. No se contempla el uso de 
este procedimiento para el diseño estruc- 
tural con bambú. El uso de las deforma- 
ciones medidas en el centro del claro de 
las piezas sometidas a flexión estática no 
genera graves errores al calcular el 
Módulo de Elasticidad del bambú, en 
lugar de hacerlo con la deformación 
máxima que se presente en otro punto 
muy cercano del centro. 
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